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ABSTRACT 

This  work  investigates  the  linear  stability  of  an  antiplane  shear  motion  which 
involves  a  planar  phase  boundary  in  an  arbitrary  element  of  a  wide  class  of  non- 
elliptic  generalized  neo-Hookean  materials  which  have  two  distinct  elliptic  phcises. 
It  is  shown,  via  a  normal  mode  analysis,  that,  in  the  absence  of  inertial  effects, 
such  a  process  is  linearly  imstable  with  respect  to  a  large  class  of  disturbances  if 
and  only  if  the  kinetic  response  function — a  constitutively  supplied  entity  which 
gives  the  normal  velocity  of  a  phase  boundary  in  terms  of  the  driving  traction 
which  acts  on  it  or  vice  versa — is  locally  decreasing  as  a  function  of  the  appro¬ 
priate  argument.  An  alternate  analysis,  in  which  the  linear  stability  problem  is 
recast  as  a  functional  equation  for  the  interface  position,  allows  the  interface  to 
be  tracked  subsequent  to  pertiu’bation.  A  particular  choice  of  the  initial  distur¬ 
bance  is  used  to  show  that,  in  the  case  of  an  unstable  response,  the  morphological 
character  of  the  phase  boundary  evolves  to  qualitatively  resemble  the  plate-like 
structures  which  are  found  in  displacive  solid-solid  phase  transformations.  In 
the  presence  of  inertial  effects  a  combination  of  normal  mode  and  energy  anal¬ 
yses  are  used  to  show  that  the  condition  which  is  necessary  and  sufficient  for 
instability  with  respect  to  the  relevant  class  of  pertmbations  in  the  absence  of 
inertia  remains  necessary  for  the  entire  class  of  perturbations  and  sufficient  for 
all  but  a  very  special,  and  physically  unrealistic,  subclass  of  these  perturbations. 
The  linear  stability  of  the  relevant  process  depends,  therefore,  entirely  upon  the 
transformation  kinetics  intrinsic  to  the  kinetic  response  function. 
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1.  INTRODUCTION 

Displacive  solid-solid  phase  transformatiors  occur  in  a  wide  variety  of  metal¬ 
lic  and  ceramic  alloys.  The  different  phases  of  a  material  capable  of  imdergoing 
such  a  transformation  are  generally  distinguished  by  distinct  crystal  structirres. 
Transformations  involving  such  materials  are  characterized  by  the  existence  of 
interfaces,  or  phase  boundaries,  which  segregate  material  in  different  phases.  Of 
particular  interest  is  the  growth  stage  of  a  displacive  solid-solid  phase  transfor¬ 
mation  which  is  directly  related  to  the  kinetics  governing  the  motion  of  phase 
boundaries.  Some  experimental  work  directed  at  understanding  the  growth  stage 
of  these  transformations  has  been  performed.  Nishiyama  [25]  has  separated  the 
transformation  kinetics  of  relevant  materialr  into  three  classes  baised  upon  the 
speed  at  which  they  occur.  Depending  on  how  they  are  loaded,  some  materials 
may  exhibit  kinetics  which  fall  into  any  or  all  of  these  three  classes.  In  the  fastest 
of  these,  phase  boundaries  propagate  at  velocities  which  are  of  the  same  order  of 
magnitude  as  the  velocity  of  shear  wave  propagation;  in  the  remaining  two  classes 
the  velocities  with  which  phase  boundaries  propagate  are  many  orders  of  magni¬ 
tude  smaller  the  the  shear  wave  speed.  The  work  of  Grujicic,  Olson  &  Gwen 
[16]  and  Clapp  &  Yu  [10]  suggests  that  slowly  propagating  phase  boundaries 
are  most  often  observed  to  be  plzuicir  in  structure,  while  those  whi''h  propagate 
rapidly  often  display  highly  complex  geometries  involving  plate- like  or  dendritic 
structures.  These  complicated  structures  are  reminiscent  of  t.iose  which  occur  in 
crystal  growth,  that  are  known  to  evolve  from  states  involving  planar  interfaces 
which  separate  solidified  crystal  material  from  liquid  melt.^  It  is,  therefore,  nat¬ 
ural  to  speculate  as  to  whether  the  complicated  plate-like  morpholgies  observed 
in  rapid  displacive  solid-solid  phaise  transformations  can  emerge  in  an  anzdogous 
fashion  from  their  slow  counterparts. 

Finite  elastic  dynamical  processes  in  materieils  capable  of  sustaining  equilib¬ 
ria  with  discontinuous  deformation  gradients  have  figured  prominently  in  recent 
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Lancer  [22]  provides  an  overview  of  such  phenomena. 
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continuum  mechanical  treatments  of  displacive  solid-solid  phase  transformations.^ 
In  a  homogeneous,  hyperelastic  material  equilibrium  states  with  discontinuous  de¬ 
formation  gradients  occur  only  if  the  relevant  elastic  potential  allows  for  a  loss 
of  ellipticity — at  certain  values  of  the  deformation  gradient — in  the  associated 
displacement  equations  of  equilibrium?  Materials  characterized  by  elastic  po¬ 
tentials  which  allow  such  a  loss  of  ellipticity  are  referred  to  as  non-elliptic.  Of 
particular  importance  in  most  of  the  work  that  has  been  done  in  this  area  are 
non-elliptic  materials  which  have  two  disjoint  elliptic  phases.  Such  hypothetical 
materials  serve  as  models  for  actual  materials  which  can  sustain  displacive  solid- 
solid  phase  transformations;  surfaces  which,  in  either  equilibrium  or  dynamics, 
separate  the  different  phases  of  a  non-elliptic  material  function  as  models  for  the 
phase  boundaries  which  occur  in  actual  materials  and,  hence,  are  referred  to  as 
such. 

Despite  the  apparent  dearth  of  experimental  information  regcirding  the  is¬ 
sue  of  whether  the  growth  stage  of  displacive  solid-solid  phase  transformations 
can  involve  the  emergence  of  complicated  dendritic  structures  from  plamar  phase 
boundaries,  it  is  legitimate  to  examine  this  topic  from  an  analytical  perspective 
in  the  foregoing  continuum  mechanical  context.  Except  for  the  work  of  Silling 
[30],  the  bulk  of  the  continuum  mechanical  investigations  which  consider  dynam¬ 
ical  processes  are  confined  to  one-dimensional  bar  theory  and,  hence,  are  not 
of  direct  bearing  on  the  issue  of  phase  boundary  morphology.  SiLLiNG  [30]  has 
demonstrated,  through  an  asymptotic  analysis,  that  a  particular  generalized  neo- 
Hookean  material  is  capable  of  sustmning  a  motion  which  involves  a  steadily  prop¬ 
agating  cusped  surface  of  discontinuity  which  segregates  distinct  elliptic  phases 
of  the  relevant  material.  This  cusped  phase  boundary  can  be  thought  of  as  a 
model  for  one  which  would  accompany  a  single  plate-like  structure  in  zm  actual 
displacive  solid-solid  phase  trzmsformation.  Silling  [30]  also  performs  numerical 


^  See,  e.g.,  Abeyaratne  &  Knowles  (4-6],  James  [18],  Pence  [27]  and  Silling  [30], 
^  For  a  discussion  of  this  issue  consult,  for  instance,  RosaKIS  [28]. 
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calculations  which  seem  to  support  his  asymptotic  results.  It  is  importeint  to  note 
that  this  work  does  not  consider  the  issue  of  the  emergence  of  the  cusped  phase 
boundary  from  a  planar  one. 

In  analogy  to  the  large  body  of  analytical  work  which  has  been  directed  at 
modeling  the  emergence  of  dendritic  structures  from  planar  interfaces  in  the  pro¬ 
cess  of  crystal  growth,'*  it  seems  reasonable — as  a  first  step  in  addressing  the  issue 
at  hand — to  investigate  the  stability  of  a  two-phase  process  involving  a  steadily 
propagating  planar  phase  boundary  in  a  non-elliptic  material.  In  a  study  which 
focuses  primarily  on  contructively  establishing  the  existence  of  two-phase  equi¬ 
libria  in  a  special  non-elliptic  generalized  neo-Hookean  material  Fried  [14]  also 
analyzes — in  an  inertia-free  setting — the  stability  of  such  a  state  with  respect  to 
a  particular  class  of  perturbations.®  Together,  the  narrow  class  of  perturbations 
which  is  considered,  the  absence  of  inertial  effects,  and  the  constitutive  special¬ 
ization  which  is  adhered  to  severely  restrict  the  generality  of  the  results  which 
are  obtained  in  [14].  The  objective  of  the  present  inquiry  is,  therefore,  to  perform 
a  more  general  stability  analysis  where  a  two-phase  process  involving  a  steadily 
propagating  planar  phase  boundary  in  a  wide  class  of  non-elliptic  generalized 
neo-Hookean  materials  is  subjected  to  a  broad  class  of  disturbances  and  inertial 
effects  are  taken  into  consideration.  It  will  transpire,  however,  that  the  stability 
results  which  eire  obtained  are  consistent  with  those  secured  by  Fried  [14]. 

Chapter  2  is  devoted  to  preliminaries.  After  a  brief  overview  of  the  notation 
to  be  used.  Section  2.1  introduces  the  kinematics  and  fundamental  balance  prin¬ 
ciples  which  will  be  needed  in  the  following.  Section  2.2  explains  the  constitutive 
restrictions  which  will  be  adhered  to  throughout  this  work.  Section  2.3  is  con¬ 
cerned  with  the  notions  of  mechanical  dissipation  and  driving  traction  which  are 
associated  with  phase  boundaries;  these  lead  naturally  to  the  consideration  of  a 
kinetic  relation — which  gives  the  normal  velocity  of  a  phase  boundary  in  terms 

■*  See,  for  example,  Langer  [22],  Mullins  &  Sekerka  [24]  and  Strain  [31], 

®  The  relevant  material  is  also  that  used  by  SiLLiNG  in  [29-30]. 
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of  the  driving  traction  which  acts  on  it  or  vice  versa — zmd  em  associated  kinetic 
response  function.  In  the  final  section  of  Chapter  2,  the  kinematics  are  specialized 
to  those  of  antiplane  shear. 

Chapter  3  concentrates  upon  a  linear  stability  analyis  of  a  two-phase  process 
involving  a  steadily  propagating  planar  phase  boimdary  in  a  non-elliptic  general¬ 
ized  neo-Hookean  material  which  obeys  the  Baker-Ericksen  inequality.  The  pro¬ 
cess  to  be  perturbed,  which  involves  an  antiplaue  shear  deformation,  is  introduced 
in  Section  3.1.  In  Section  3.2  the  class  of  perturbations  which  will  be  applied  to 
the  base  process  are  then  introduced.  Each  admissible  perturbation  involves,  in 
general,  a  disturbance  of  the  configuration  of  the  phase  boundary  and  of  the  dis¬ 
placement  and  velocity  fields  in  a  small  neighborhood  of  the  phase  boundary — all 
of  which  eire  assumed  to  be  small  in  some  appropriate  sense.  The  kinematics  of 
the  perturbation  axe  also  restricted  to  those  of  antiplane  shear.  It  is  assumed, 
furthermore,  that  the  post-pertiurbation  deformation  remains  an  antiplane  shear 
and  involves  only  one  phase  boundary.  Section  3.3  is  devoted  to  the  linearization 
about  the  base  process  of  the  field  equations,  which  hold  away  from  the  phase 
boundary,  about  the  base  process  introduced  in  Section  3.1.  In  a  similar  manner, 
Section  3.4  is  concerned  with  the  lineeirization  about  the  beise  process  of  the  jiunp 
conditions  and  kinetic  relation  which  hold  on  the  phase  boundary.  A  summary 
of  the  complete  linearized  system  of  field  equations,  jump  conditions,  kinetic  re¬ 
lation,  boundary  and  initial  conditions  which  describe  the  process  generated  by 
the  perturbation  is  presented  in  Section  3.5.  Included  are  both  the  inertial  eind 
inertia-free  czises.  In  Section  3.6  a  normal  mode  analysis  is  performed  in  the  ab¬ 
sence  of  inertial  effects.  A  condition  necessaxy  and  sufficient  for  the  base  process 
to  be  unstable  with  respect  to  any  perturbation  of  the  type  introduced  in  Section 
3.2  is  obtzuned.  This  condition  involves  only  the  local  behavior  of  the  derivative  of 
the  (essentially  arbitrzuy)  kinetic  response  function  introduced  in  Section  2.3.  An 
alternative  to  the  normeil  mode  analysis  of  Section  3.6  is  performed  in  Section  3.7. 
Here  the  relevant  initial  boimdaxy  value  problem  is  converted  into  a  functional 
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initial  value  problem  for  the  correction  to  the  interface  position  arising  from  the 
perturbation.  Analysis  of  this  problem  yields  identical  stability  criteria  to  those 
which  are  obtained  in  Section  3.6.  Moreover,  when  instability  is  present,  its  man¬ 
ifestation  can  be  tracked  over  a  finite  time  interv2il  upon  which  the  linearization 
performed  in  the  foregoing  remains  valid.  The  results  suggest  the  emergence  of 
plate-like  or  dendritic  structures.  Section  3.8  contains  both  a  normal  mode  and 
an  energy  analyis  in  the  case  where  inertial  effects  are  included.  First,  the  normal 
mode  analysis  leads  to  a  necessary  and  sufficient  condition  for  the  base  process 
to  be  unstable  with  respect  to  a  special  subset  of  the  class  of  perturbations  intro¬ 
duced  in  Section  3.2.  The  operative  stability  criterion  is  identical  to  that  which 
holds  for  all  initial  disturbances  in  the  inertia-free  case.  An  argument  based  upon 
a  Fourier-Laplace  transform  analysis  of  the  relevant  initicil  boimdary  value  prob¬ 
lem  is  then  used  to  show  that  all  but  a  very  special  subset  of  the  full  class  of 
perturbations  introduced  in  Section  3.2  are  covered  by  the  normal  mode  analysis. 
Section  3.8  is  concluded  with  an  energy  analysis  which  is  used  to  show  that  the 
sufficient  condition  referred  to  above  is  also  necessary  for  the  base  process  to  be 
unstable  with  respect  to  any  perturbation  of  the  cleiss  introduced  in  Section  3.2. 
Hence,  it  is  shown  that  the  presence  of  inertia  does  not  qualitatively  alter  the 
linear  stability  of  the  base  process  of  interest.  Finally,  in  Section  3.9,  a  discus¬ 
sion  which  focusses  on  the  physical  reasonableness  of  admissible  non-monotonic 
kinetic  response  functions  is  undertaken. 
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2.  PRELIMINARIES 

2.1.  Notation,  kinematics  and  balance  principles.  In  the  following  IR 

and  (T  denote  the  sets  of  real  and  complex  numbers.  The  intervals  (0,  oo)  and 

[0,  oo)  are  represented  by  j?2+  and  M+.  The  symbol  M",  with  n  equal  to  2  or  3, 

represents  real  n-dimensional  space  equipped  with  the  standard  Euclidean  norm. 

—  ® 

If  17  is  a  set,  then  its  closure,  interior  and  boundary  are  designated  by  Z7,  U ,  and 
dU,  respectively.  The  complement  of  a  set  V  with  respect  to  U  is  written  as 
U  \  V.  Given  a  function  ip  :  U  ^  W  and  a  subset  V  of  U,  ip{V)  stands  for  the 
image  of  V  under  the  map  ip. 

Vectors  and  linear  transformations  from  to  ]R^  (referred  to  herein  as 
tensors)  are  distinguished  from  scalars  with  the  aid  of  boldface  type — lower  and 
upper  case  for  vectors  and  tensors,  respectively.  Let  a  and  b  be  vectors  in  iR^, 
their  inner  product  is  then  written  as  a  •  b;  the  Euclidean  norm  of  a  is,  further, 
written  as  |a|  =  y/a.  •  a.  The  set  of  unit  vectors — that  is,  vectors  with  unit 
Euclidean  norm — in  ]R^  is  designated  by  M.  The  symbol  C  refers  to  the  set 
of  tensors,  £+  denotes  the  set  of  all  tensors  with  positive  determinant,  and  S 
stands  for  the  collection  of  all  symmetric  positive  definite  tensors.  If  F  is  in  £ 
then  represents  its  transpose;  if,  moreover,  det  F  ^  0,  then  the  inverse  of  F 
and  its  transpose  axe  written  as  F~*  and  F“^,  respectively.  The  notation  a®  b 
refers  to  the  tensor  A,  formed  by  the  outer  product  of  a  with  b,  defined  such  that 
Ac  =  (b  •  c)a  for  any  vector  c  in  iR^.  If  A  and  B  are  tensors  then  their  inner 
product  is  written  as  A  •  B  =  tr  AB^. 

When  component  notation  is  used,  Greek  indices  range  only  over  {1,2}; 
summation  of  repeated  indices  over  the  appropriate  range  is  implicit.  A  subscript 
preceded  by  a  comma  denotes  partial  differentiation  with  respect  to  the  corre¬ 
sponding  coordinate.  Also,  a  superposed  dot  signifies  partial  differentiation  with 
respect  to  time. 

Consider,  now,  a  body  B  which,  in  a  reference  configuration,  occupies  a  region 
Tl.  contained  in  IR^.  A  motion  of  6  on  a  time  interval  T  C  iR  is  characterized  by 
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a  one-parameter  family  of  invertible  mappings  y{-,t)  :  TZ  —*  TZt,  with 

y(x,<)  =  X -f  u(x,<)  V(x,t)GX,  (2.1.1) 

where  Ad  —  TZ  x  T  represents  the  trajectory  of  the  motion.  Assume  that  the 
deformation  y,  or  equivalently  the  displacement  u,  is  continuous  and  possesses 
piecewise  continuous  first  and  second  partial  derivatives  on  A4.  Let  St  be  the  set 
of  points  contained  in  TZ  defined  so  that,  at  each  instant  t  in  T,  y(-,t)  is  twice 
continuously  differentiable  on  the  set  'JZ\  St.  Let  the  set  E  be  defined  by 

E  =  {(y.,t)\xeSt,teT}.  (2.1.2) 

Introduce  the  deformation  gradient  tensor  F  :  Ai  \  E  C  hy 

F(x,  t)  =  Vy(x,  t)  V(x,  t)  6  X  \  r,  (2.1.3) 

where  the  associated  Jacobian  determinant,  J  :  M  \  E  M,  of  y  is  restricted  to 
be  strictly  positive  on  its  domeiin  of  definition: 

7(x,t)  =  detF(x,t)  >  0  V(x,<)  6  Ad  \  17. 

+ 

Hence,  F  :  Ad  \  i7  — >  C+.  The  left  Cauchy- Green  tensor  G  :  Ad  \  Z"  — >  <5 
corresponding  to  the  deformation  y  is  given  by 

G(x,  t)  =  F(x,  t)F^(x,  t)  V(x,  t)eM\E.  (2.1.4) 

The  deformation  invariants  associated  with  y  exist  on  Ad  \  Z  and  Eire  supplied 
through  the  fundamented  scalar  invariants  of  G; 

7i(G)  =  trG,  J2(G)  =  i((trG)2-tr(G2)),  /3(G)  =  detG.  (2.1.5) 


-s- 

With  the  above  kinematic  antecedents  in  place  introduce  the  nominal  mass 
density  p  the  nominal  body  force  per  unit  mass  b  :  and  the 

nominal  stress  tensor  S  :  Ad  \  X*  — >  £,  and  suppose  that  p  is  constant  on  %  and  b 
is  continuous  on  Ad,  while  S  is  piecewise  continuous  on  Ad,  continuous  on  Ad  \  X, 
and  has  a  piecewise  continuous  gradient  on  Ad.  Let  p*  be  the  mass  density  in 
the  deformed  configiuration  associated  with  y.  Given  a  regular  subregion  V  of 
71,  with  87^  n  5f  a  set  of  measure  zero  in  67^  for  each  <  in  T,  let  m  :  dV  M 
denote  the  unit  outward  normal  to  dV.  Then  the  global  balance  laws  of  mass, 
linear  momentum,  and  angular  momentum  require  that 


on  T, 


(2.1.6) 


respectively,  for  every  such  regular  subregion  7^  contained  in  7Z. 

Localization  of  the  balance  laws  (2.1.6)-(2.1.8)  at  an  arbitraiy  point  con¬ 
tained  in  the  interior  of  Ad  \  X  yields  the  following  familiar  field  equations: 


P  =  P*{y)J  on  Ad  \  X, 

V  •  S -f  pb  =  pii  on  Ad  \  X,  (2.1.9) 

SF^=FS'^  on  Ad\X. 


Suppose,  from  now  on,  that  the  set  St  is  a  reguleir  surface  for  every  t  in  T. 
The  set  X  then  represents  the  trajectory  of  a  surface  of  discontinuity  in  F  and  S. 
Let  g{-,t)  denote  a  generic  field  quantity  g{-,t)  :  St  IR  which  is  discontinuous 
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across  St  at  the  insteint  t  in  T.  Define  the  jump  of  across  St  by 


l9(x,t)I  =  lim  (sf(x  +  hn(x,t),<) -gr(x  - /in(x, <),<))  V(x,t)€r,  (2.1.10) 


where  n(-,t)  :  5t  —»  A/"  is  a  unit  normal  to  St  at  each  t  in  T.  Then,  localization 
of  (2.1.6)-(2.1.8)  at  an  arbitrary  point  in  U  yields  the  following  jump  conditions 


[p,(y)J|  =  0  on  E, 
[Sn]  +  pV„lu|  =  0  on  E, 


(2.1.11) 


where  Vn{-,t)  :  St  JR  denotes  the  component  of  the  velocity  of  the  surface  St 
in  the  direction  of  n(^-,t)  at  the  instant  t  in  T. 

Equations  (2.1.9)i  and  (2.1.11)i  are,  evidently,  completely  decoupled  from 
equations  (2. 1.9)2, 3  and  (2.1.11)2;  that  is,  given  a  solution  to,  say,  a  boundary 
value  problem  involving  (2. 1.9)2, 3  and  (2.1.10)2,  p*  can  be  cedculated  a  posteri- 
ori.  For  this  reason  equations  (2.1.9)i  and  (2.1.11)i  will  be  disregarded  in  the 
subsequent  analysis. 

In  this  investigation  an  inertia-free  motion  is  defined  as  one  wherein  the 
inertial  terms  on  the  right  hand  sides  of  the  global  balance  equations  (2.1.7)  and 
(2.1.8)  are  replaced  by  the  zero  vector.  In  the  context  of  an  inertia-free  motion 
the  field  equation  (2. 1.9)2  simplifies  to  read 


V  •  S  -t-  pb  =  0  on  M  \  E, 


(2.1.12) 


and  the  jump  condition  (2.1.11)2  becomes 

|Snl=0  on  E.  (2.1.13) 

Equations  (2.1. 9)1,3  and  (2.1.11)i  remain,  of  course,  unaltered. 
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In  addition  to  the  jump  conditions  given  in  (2.1.11)  in  the  inertial  case  or 
(2.1.10)1  and  (2.1.13)  in  the  inertia-free  case,  the  stipxilated  continuity  of  y  gives 
the  following  kinematic  jump  condition 

fu|  =0  on  r.  (2.1.14) 

2.2.  Constitutive  assumptions.  Let  B  be  composed  of  a  hyperelastic  ma¬ 
terial  which  is  homogeneous,  isotropic  and  incompressible.  Since  B  is  hyperelastic 
its  mechanical  response  is  governed  by  zin  elastic  potential  or  strain  energy  per 
unit  reference  volume.  The  homogeneity  of  B  implies  that  the  elastic  potential 
does  not  depend  explicitly  on  position  in  the  reference  configuration.  Further¬ 
more,  because  B  is  isotropic  the  elastic  potential  can  depend  on  the  deformation 
gradient  F  only  through  the  deformation  invariants  Ik(G)  defined  in  (2.1.5).  The 
incompressibility  of  B  requires  that  the  deformation  y  be  isochoric.,  t.e., 

/3  (G(x,  t))  =  J^(x,  t)  =  l  V(x,  <)  €  M  \  r.  (2.2.1 ) 

An  additional  consequence  of  isotropy  is,  therefore,  that  the  elastic  potential  can 
be  expressed  as  a  function  solely  of  the  first  two  deformation  invariants.  It  can 
also  be  demonstrated  via  (2.1.5)  that,  when  (2.2.1)  holds,  /q(G(x,<))  >  3  for  all 
(x,  f)  contained  in  A4  \  X".  Now,  let  W  :  [3, 00)  x  [3, 00)  JR  denote  an  elastic 
potential  which  characterizes  B  and  zissume  that  W  is  continuously  differentiable 
with  piecewise  continuous  second  derivatives  on  its  domain  of  definition.  The 
nominal  stress  response  of  B  is  then  determined  through  IF  up  to  an  arbitrary 
pressure  p  :  M\I!  JR  required  to  accomodate  the  kinematic  constraint  (2.2.1) 
imposed  by  the  incompressibility  of  B:  viz., 

S  =  2  (lF/,(/)F  -h  Wi,{I){li{G)l  -  G)f)  -pF-^  on  Af  \  T,  (2.2.2) 
where  I :  M\S  ^  13,oo)  x  (3,oo)  is  given  by 


/(x,()=  (/,(G(x,()),/2(G(x,()))  V(x,()6 
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Following  Gurtin  [17],  let  the  class  of  generalized  neo-Hookean  materials  re¬ 
fer  to  that  subset  of  hyperelastic  materials,  first  introduced  by  Knowles  [19], 
which  are  homogeneous,  isotropic  and  incompressible  with  elastic  potential  in¬ 
dependent  of  the  second  deformation  invariant  (2. 1.5)2.  Assume,  henceforth, 
that  B  is  composed  of  a  generalized  neo-Hookeain  material  with  elastic  potential 
W  :  [3,  oo)  — »  jR,  where  W  is  continuously  differentiable  with  piecewise  contin¬ 
uous  derivative  on  [3,  oo).  Then,  by  (2.2.2),  the  nominal  stress  response  of  B  is 
determined  by 

S  =  21F'(/i(G))F-pF-^  on  M\E.  (2.2.3) 

Suppose  also  that  the  elastic  potential  is  normalized  so  that 

W^(3)  =  0.  (2.2.4) 

Choose  a  rectangular  Cartesian  frame  X  =  {0; 61,62, 63}  and  consider  the 
response  of  the  material  at  hand  to  a  simple  shear  deformation  y  given  by 

y(x,<)  =  (14- 763  0  ei)x  V(x,t)€Af,  (2.2.5) 

where  the  constant  7 — assumed  non-negative  without  loss  of  generality — denotes 
the  amount  of  shear.  From  (2.1.3),  (2.2.3)  and  (2.2.5)  the  nominal  shear  stress 
corresponding  to  the  deformation  y  is,  for  each  7  in  iR+,  found  to  be 

63  •  S6i  =  2jW'i3  +  7^)  =:  r(7).  (2.2.6) 

In  [19-20]  Knowles  demonstrates  that  the  31  and  32  components  of  nom¬ 
inal  and  Cauchy  shear  stress  are,  in  the  present  setting,  equal.  The  function 
r  :  1R+  —*■  IR  is,  hence,  referred  to  as  the  shear  stress  response  function  of  the 
generalized  neo-Hookean  material,  characterized  by  W,  in  simple  shear.  An  im¬ 
mediate  consequence  of  (2.2.4)  and  (2.2.6)  is 

V  /i— 3 

W{h)=  J  T{j)dj  V/iE[3,oo), 

0 


(2.2.7) 
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so  that  the  response  of  a  generalized  neo-Hookean  matericil,  in  all  three  dimen¬ 
sional  deformations,  is,  up  to  a  hydrostatic  pressure,  completely  characterized 
by  specifying  the  shear  stress  response  function  r.  Define  the  secant  modulus  in 
shear  M  :  jR+  — »  jR  of  a  generalized  neo-Hookean  material  with  elastic  potential 
W  by 

M(7)  =  2T^'(3  -1-7^)  V7  6  ^+,  (2.2.8) 

and  assume  that,  in  compliance  with  the  Baker- Ericksen  inequality, 

M(7)>0  VyeM+.  (2.2.9) 

Assume,  also,  that  M(0)  >  0  so  that  the  infinitesimal  shear  modulus  of  the 
material  at  hand  is  positive.  Note  from  (2.2.6)  and  (2.2.8)  that  the  shear  stress 
response  function  r  must  also  satisfy 

r(0)  =  0,  r'(0)  =  M(0).  (2.2.10) 

Observe,  in  addition,  that  the  stipulated  smoothness  of  W  guarantees  that  both 
r  and  M  are  piecewise  continuously  differentiable  on  iR+. 

Despite  the  significant  restrictions  which  have  been  placed  upon  the  class  of 
materials  which  will  be  considered  in  this  investigation,  no  presuppositions  have 
been  made  regarding  the  sign  of  the  derivative — where  it  exists — of  the  shear 
stress  response  function  corresponding  to  the  generalized  neo-Hookean  material 
defined  through  (2.2.6).  In  [20]  Knowles  shows  that  the  monotonicity  of  the 
shear  stress  response  function  r  is  related  directly  to  the  ellipticity  of  the  gener¬ 
alized  neo-Hookean  material  which  it  characterizes:  if  r  is  not  a  monotonically 
increasing  fimction  on  its  domain  of  definition  then  the  eissociated  material  is  non- 
elliptic.  This  investigation  will  make  use  of  a  particular  subclass  of  non-elliptic 
generjilized  neo-Hookean  materials,  first  suggested  by  Abeyaratne  [l];  this  class 
of  materials  is  characterized  by  the  set  of  shear  stress  response  functions  r  which 


are  continuous  on  M+  and  piecewise  continuously  differentiable  on  \  {7,7}, 

where  0  <  7  <  7,  such  that 

* 


t'>0  on  iR+\[7,7], 

* 

t'  <  0  on  (7,7). 


(2.2.11) 


The  sets  of  shear  strains  lying  in  the  intervals  [0,7)  and  (7,  00)  axe  referred 
to  as  the  high  and  low  strain  phases  of  the  generalized  neo-Hookeain  material 
specified  by  the  shear  stress  response  function  r.  Together  the  high  and  low 
strain  phases  of  such  a  material  comprise  its  elliptic  phases.  A  generalized  neo- 
Hookean  material  characterized  by  a  shear  stress  response  function  of  this  type 
will  be  referred  to  herein  as  a  three-phase  material.  See  Figure  1  for  a  graph  of  a 
shear  stress  response  function  typical  of  those  which  specify  three-phcise  materials. 


'  2.3.  Dissipation,  driving  traction  and  the  kinetic  relation.  Let  V  be 
a  regular  subregion  contained  in  TZ  chosen  so  that  dV  H 17  is  a  set  of  measure  zero 
in  dV.  The  total  mechanical  energy  contained  in  P  at  an  instant  t  contained  in 
T  is  given,  under  the  present  constitutive  assumptions,  by 


E(t;V)  =  +  'it&T.  (2.3.1) 

V 

A  st2mdajd  calculations  then  shows,  with  the  aid  of  (2.2.3),  that 

E{t-,'P)  =  /  S(x,<)n(x,<)  •  u(x,t)dA  -  j  /(x,<)V'„(x,f)dA  ^teT,  (2.3.2) 
dv  VnS, 

where  f  :  E  —*  JR  is  the  scalar  driving  traction  given  by 


/(•,()  =  ItV(/.(G(.,()))l-((S(-,()))-(F(.,()l  on  5.  V(€T,  (2.3.3) 


and,  given  a  generic  field  quantity  g(-,t)  :  St  —*  JR  which  jumps  across  St  at  the 
instant  t  in  T,  ((fl'(-,<)))  is  defined  through 
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In  [21]  Knowles  derives  (2.3.2)  in  the  absence  of  inertial  effects;  following  either 
Yatomi  k  Nishimura  [32]  or  Abeyaratne  k  Knowles  [3]  it  can  be  shown  that 
(2.3.3)  reduces,  in  this  setting,  to 

/(•,t)  =  (W(Ji(G(-,t)))l-S(-,<)-lF(-,01  on  5*  VteT,  (2.3.5) 

where  S(-,f)  {resp.,  S(-,t))  is  the  limiting  value  of  the  field  S(-,<)  on  the  side  of 
the  interface  into  which  the  vmit  normal  n(-,<)  is  (resp.,  is  not)  directed  at  the 
instant  t  in  T. 

From  (2.3.2)  it  is  clear  that  the  presence  of  a  moving  surface  of  discontinuity 
St  of  the  type  considered  here  may  effect  the  balance  of  mechanical  energy.  Let  the 
difference  in  the  rate  of  work  of  the  mechanical  forces  external  to  V  and  the  rate  at 
which  energy  is  stored  in  V  be  referred  to  cis  the  rate  of  dissipation  of  mechanical 
energy  associated  with  the  region  V.  When  treated  from  a  thermomechanical 
perspective,  the  dissipation  rate  can  be  shown  to  be  identical  to  the  product  of  the 
temperature  and  the  rate  of  entropy  production — provided  that  the  temperature  is 
spatially  uniform  and  independent  of  time.®  The  Clausius-Duhem  inequality  then 
requires  that  the  dissipation  rate  associated  with  a  motion  of  the  kind  envisioned 
here  be  non-negative,  t.e., 

I  f(y.,t)Vnix,t)dA>0  yteT,  (2.3.6) 

VnS, 

for  every  regular  subregion  V,  with  dV  H  5t  a  set  of  measure  zero  in  dV  at  t  in 
T,  contained  in  71.  A  localization  of  (2.3.6)  at  an  arbitrary  point  on  the  interface 
therefore  yields  the  inequality 

fVn  >  0  on  r  (2.3.7) 

as  a  condition  imposed  for  the  admissibility  of  the  motion. 

®  For  a  detailed  discussion  of  the  these  issues  see  Abeyaratne  &  Knowles  [3]. 
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Assume  that  the  surface  St  separates,  for  each  t  in  T,  low  and  high  strain 
elliptic  phases  of  the  three  phase  material  at  hand.  In  the  context  of  a  motion 
which  involves  such  an  interface  it  is  necessary  (see  [3-6])  to  supplement,  in  some 
fashion,  the  constitutive  information  which  relates  the  stress  and  strain  fields.  An 
approach  to  this  taken  by  Abeyaratne  &  Knowles  [3]  entails  the  provision  of 
a  kinetic  relation  which  gives  the  normal  velocity  of  the  interface  in  terms  of  the 
driving  traction  that  acts  on  it  or  vice  versa.  In  the  former  case  a  constitutive 
fimction  V  :  M  —*  IR  is  provided  so  that 

Vn  =  V{f)  V/  e  M,  (2.3.8) 

while,  in  the  latter  case  a  constitutive  function  :  M  —*  M  is  given  so  that 

/  =  <^(V„)  VK  €  M.  (2.3.9) 

The  functions  V  and  are  referred  to  as  the  kinetic  response  functions.  Both 
varieties  of  kinetic  response  functions  will  be  considered  in  this  investigation.  If 
V  is  such  that  V{f)f  >  0  for  cdl  /  in  iR  then  (2.3.8)  is  automatically  satisfied 
and  V  is  referred  to  as  admissible.  If  (f>{V)V  >  0  for  all  V  in  JR,  ip  is,  similarly, 
referred  to  as  admissible.  If  an  admissible  kinetic  response  function  V  (or  ip)  is 
continuous  on  JR,  then  it  must  satisfy  F(0)  =  0  (or  (,5(0)  =  0).  If,  in  addition,  to 
being  admissible,  V  (or  (f>)  is  continuously  differentiable  on  JR,  then  K'(0)  >  0  (or 
(^'(0)  >  0).  Otherwise,  admissiblility  implies  nothing  with  regard  to  the  sign  of  the 
derivative  of  a  smooth  kinetic  response  function.  All  kinetic  response  functions 
considered  herein  are  assumed  to  be  admissible.  See  Figure  2  and  Figure  3  for 
graphs  of  such  kinetic  response  functions. 

Abeyaratne  [1],  Ball  &  James  [9],  Gurtin  [17],  Fosdick  h  MacSithigh 
[12],  and  Silling  [29]  consider  either  equilibrium  states  or  inertia-free  motions 
and  require  that  the  driving  traction,  /,  be  identically  equal  to  zero  on  E.  This  is 
equivalent  to  prescribing  a  supplementary  kinetic  relation  in  the  form  (2.3.9)  with 


If  identically  zero  on  IR.  Provision  of  such  a  kinetic  relation  is,  furthermore,  a 
necessary  consequence  of  requiring  that  a  suitable  energy  functional  be  minimized 
at  each  <  in  T  (see  Abeyaratne  [2]). 

2.4.  Antiplane  shear  motions  of  a  generalized  neo-Hookean  mate¬ 
rial.  Suppose,  from  now  on,  that  TZ  is  &  cylindrical  region  and  choose  a  rectangu¬ 
lar  Cartesian  frame  X  =  {0; €1,02,63}  so  that  the  imit  base  vector  63  is  parallel 
to  the  generatrix  of  TZ.  The  deformation  y  defined  through  (2.1.1)  consists  of  an 
antiplane  shear  normed  to  the  plane  spanned  by  the  base  vectors  ei  and  02  if  it 
is  of  the  form 

y{-K,t)  =  X  +  u{xi,x2,i)e3  V(x,<)  6  Af.  (2.4.1) 

Observe  that  the  displacement  field  intrinsic  to  an  antiplane  shear  deformation 
of  this  type  has  only  one  nonzero  component  which  lies  in  the  03  direction  and 
is  independent  of  the  13-coordinate.  In  (2.4.1)  Xq  =  x  •  eo  for  each  x  contained 
in  TZ.  The  function  u  will  be  referred  to  as  the  out-of-plane  displacement  field. 
Inspection  of  (2.4.1)  reveals  that  any  discontinuities  in  the  gradient  of  y  must 
result  from  discontinuities  in  the  out-of-plane  displacement  field  and,  hence,  occur 
across  surfaces  which  do  not  vary  with  the  X3-coordinate.  Let  St  denote  such  a 
surface  at  the  instant  t  in  T  emd  let  X  be  defined  as  in  (2.1.2). 

It  is  possible  to  show,  following  the  work  of  Knowles  [20]  in  the  inertia-free 
context,  that,  edthough  not  every  hyperelastic,  isotropic  and  incompressible  ma- 
tericd  can  sustain  nontrivial  antiplane  shear  motions,  all  generalized  neo-Hookean 
materials  are  capable  of  doing  so.  It  is  easily  shown  that  for  such  materials  the  lo- 
ced  balance  equations  (2. 1.9)2, 3  reduce,  in  the  absense  of  body  forces  and  under  the 
asssumption  that  the  nominal  stress  tensor  is  independent  of  the  X3 -coordinate, 
to  the  scalar  equation 


{M{j)u,a),a  =  pii  on  X\r,  (2.4.2) 

where  X  is  given  hy  V  y.  T,  V  is  a  generic  cross  section  of  TZ,  and  F  = 
{(ii,X2,<)l(xi,X2)  €  Ct,t  E  T}  with  Ct  =  P  n  5(  at  each  t  in  T.  See  Fosdick 
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&  Serrin  [13]  zind  Fosdick  &  Kao  [11]  for  a  general  discussion  of  circumstances 
under  which  the  field  equations  (2. 1.9)2, 3  reduce  to  a  single  scalar  equation.  In 
(2.4.2)  M  is  the  secant  modulus  in  shear  as  defined  in  (2.2.8)  eind  7  :  A"  \  jT  1R 
is  the  shear  strain  field  given  by 

7(xi,i2,t)  =  V u,a{xi,X2,t)u,c,ixi,X2,t)  £  X  \  F.  (2.4.3) 

For  a  generalized  neo-Hookean  material  subjected  to  antiplane  shear,  the 
jump  condition  (2.1.11)2  reduces  to 


[M(7)u,q  Haj  +  pV„[u]  =  0  on  r, 


(2.4.4) 


where  F  =  {(x,t)|x  €  Ct^i  €  T}  and  n(-,f)  :  Ct  — v  is  a  unit  normal  to  Ct, 

while  the  kinematic  jump  condition  (2.1.14)  becomes 


(u]=0  on  r. 


(2.4.5) 


It  is  also  readily  shown  that  the  driving  traction  /,  introduced  in  Section  2.3, 
for  a  generalized  neo-Hookean  material  subjected  to  an  antiplane  shear  deforma¬ 
tion  involving  a  discontinuity  in  the  gradient  and,  perhaps,  the  peirtial  derivative 
with  respect  to  time  of  out-of-plane  displacement  field  across  a  moving  curve  Ct 
is  given,  with  the  aid  of  (2.3.5),  by 


=  J_  ’■(7)<iT-  ((JW(7)u,o  ))[“,.  1 


(2.4.6) 


With  reference  to  (2.1.12),  (2.1.13)  and  (2.3.5)  it  is  easily  demonstrated  that, 
in  the  absence  of  inertial  effects,  (2.4.2)  is  replaced  by 


(M(7)u,o),a  =  0  on  A’\r, 


(2.4.7) 
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(2.4.8) 

f  =  J_  T(j)dy  -  M{y)u,c,[u,a}  on  F.  (2.4.9) 

Observe  that,  within  the  context  of  an  antiplane  shear  deformation  of  the 
type  described  above,  no  generality  is  lost  by  focussing  exclusively  upon  the 
motion  on  a  cross-section  D  of  the  cylinder  TZ  and  the  dynamics  of  the  curve 
Ct  =  F>  r\  St-  In  the  following,  curves  Ct  across  which  the  gradient  and,  perhaps, 
the  partial  derivative  with  respect  to  time  of  the  out-of-plane  displacement  field 
u(-,  •,  t)  jump,  at  some  instant  t  in  T,  and  which  segregate  the  high  and  low  strain 
phases  of  the  material  at  hand  will,  therefore,  be  referred  to  as  phase  boundaries. 


while  (2.4.4)  becomes 

[M(7)u,a  UcJ  =  0  on  r, 

and  (2.4.6)  simplifies  to 
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3.  LINEAR  STABILITY  OF  A  PROCESS  INVOLVING  A  STEADILY 
MOVING  PLANAR  PHASE  BOUNDARY  IN  A  THREE-PHASE 
MATERIAL 


3.1.  Description  of  the  base  process.  Suppose  that  B  is  composed  of 
a  three-phase  material  and  that  the  cylinder  72.  degenerates  so  as  to  occupy  all 
of  IR^.  Let  the  rectangular  Cartesian  frame  X  be  as  in  Section  2.4.  Consider 
an  antiplane  shear  motion  on  the  time  interval  (— oo,0)  with  an  out-of-plane 
displacement  field  uo(-,i)  :  IR  IR  given  by 


7jii  -f  vit 
7rri  +  Vrt 


if 

if 


a:i  <  uot, 

xi  >  Vat, 


(3.1.1) 


for  each  i  in  (— oo,  0),  where  the  shear  strains  7/  and  7r  satisfy  one  of  the  following: 


0  <  7r  <  7  <  7  <  7/,  0  <  7/  <  7  <  7  <  7r.  (3.1.2) 

*  « 

Since  one  of  (3.1.2)  must  hold,  there  is  no  loss  in  generality  incurred  by  assuming 
that  the  base  interface  normal  velocity  vq  is  non-negative;  that  is, 


vq  >  0.  (3.1.3) 

It  is  clear  that  uq  satisfies  the  differential  equation  in  (2.4.2)  on  the  set 
{IR^  X  (-oo,0))\ro  with  To  given  by  {(n,  12,  <)l  (^112:2)  6  At,<  6  (-00, 0)}  and 
At  =  {(xi,  12)1  3:1  =  vot,X2  €  IR}  for  each  t  in  (— oo,0).  The  moving  line  At  is, 
for  each  t  in  (—00, 0),  a  phase  boundary. 

Assume,  in  order  to  comply  with  the  jump  conditions  in  (2.4.4)  and  (2.4.5)  on 
Fo,  that  the  constants  7/,  7r,  vi,  Vr,  and  vo  associated  with  (3.1.1)  are  restricted 
to  satisfy  the  following  equations; 


Vr  -  VI  +  Voijr  -  Jl)  =  0, 
T(7r)  -  +  pVo{Vr  -  Vl)  =  0. 


(3.1.4) 
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Assume  that  the  normal  velocity  of  the  phase  boundary  in  the  base  process 
is  locally  subsonic  so  that  uq  satisfies  the  following  inequality: 

uo  <  min  (3.1.5) 

It  is  then  permissible^  to  impose  the  kinetic  relation  of  the  form  (2.3.8)  or  (2.3.9) 
on  Fo  and  require  that  the  parameters  7/,  7r,  u/,  tv,  and  vo  satisfy  one  of 


1^0  =  V{fo),  fo  =  ‘pM,  (3.1.6) 

depending  upon  whether  a  kinetic  relation  of  the  form  (2.3.8)  or  (2.3.9)  is  pro¬ 
vided.  In  (3.1.6)  the  base  driving  traction  fo  is  given,  with  the  aid  of  (2.4.6), 
by 

tr 

fo  =  J  r(7)  d-y  -  •|(r(7r)  +  r(7/))(7r  -  7/).  (3.1.7) 

7/ 

Observe,  as  a  consequence  of  (3.1.3)  and  (2.3.7),  that  fo  must  satisfy 

/o  >  0.  (3.1.8) 

In  a  coordinate  frame  moving  with  the  phase  boundary,  the  base  process 
described  involves  a  piecewise  homogeneous  shear  strain  field.  If  7/  and  7r  are 
consistent  with  (3.1.2)i  then  (3.1.3)  implies  that  the  base  process  is  one  wherein 

the  high  strain  elliptic  pheise  of  the  material  at  hand  grows  at  the  expense  of  the 

low  strain  elliptic  phase;  whereas,  if  -yi  and  7r  comply  with  (3. 1.2)2  then  (3.1.3) 
imphes  that  the  base  process  is  such  that  the  low  strain  elliptic  phase  of  the 
material  at  hand  grows  at  the  expense  of  the  high  strain  elliptic  phase.  In  either 
case  the  discontinuity  involved  is,  for  the  duration  of  the  motion,  a  normal  phase 
boundary — that  is,  the  angle  between  the  limiting  values  of  the  gradient  of  the 

^  See  Abeyaratne  &  Knowles  [4]. 
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out-of-plane  displacement  field  on  either  side  of  the  phase  boimdary  is  zero  at 
every  point  of  the  phase  boundary  over  the  time  interval  (— oo,0). 

Suppose,  in  addition  to  all  the  above,  that  the  kinetic  response  function  V 
or  (f  is  chosen  so  that  its  derivative  is  non-zero  at  the  base  driving  traction  /o; 
that  is,  assiune  that  one  of  the  following — as  is  appropriate  to  either  (2.3.8)  or 
(2.3.9) — must  hold  : 

V'(fo)  ^  0,  <^'(uo)  #  0.  (3.1.9) 

This  assvunption  is  made  in  order  to  preclude  the  necessity  of  going  to  higher 
order  in  the  context  of  the  forthcoming  stability  analysis.  See  Figure  2  for  the 
graph  of  a  smooth  admissible  kinetic  reponse  function  which  satisfies  (3.1.9). 

When  inertial  effects  are  ignored  it  is  clear  that  uo  as  defined  in  (3.1.1)  also 
satisfies  the  field  equation  in  (2.4.7)  on  (IR^  x  (— oo,0))  \  To-  Equation  (3.1.4)i 
is,  in  this  context,  still  sufficient  to  satisfy  (2.4.5)  on  Fq.  In  place  of  (3. 1.4)2,  the 
constants  7/,  7r,  vj,  Vr,  and  vq  must,  however,  satisfy 

^(7r)  -  t(7j)  =  0,  (3.1.10) 

in  order  for  the  jump  condition  in  (2.4.8)  to  hold  on  Jo-  Although  the  expression 
for  the  base  driving  traction  fo  given  in  (3.1.7)  remains  valid  in  the  inertia-free 
setting,  (3.1.10)  can  be  used  to  show  that 

7r 

/o  =  yr(7)d7-T,(7r-7/),  (3.1.11) 

7l 


where  r,  =  r(7/)  =  T(7r). 

Given  a  sheeir  stress  response  function  t  which  describes  a  particular  three- 
phase  material  and  an  arbitrary  kinetic  response  fimction  V  ov  (p  which  describes 
the  dynamics  of  phase  boundaries  which  may  occur  therein,  there  may  or  may 
not,  in  general,  exist  constants  7/,  7r,  uj,  Vr,  and  vq  which  satisfy  one  of  (3.1.2)i 
or  (3. 1.2)2  and  zure  consistent  with  the  restrictions  embodied  by  (3.1.4),  (3.1.5), 
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(3.1.6)i  or  (3. 1.6)2,  and  (3.1.9)i  or  (3. 1.9)2,  or  in  the  inertia-free  ccuse,  (3.1.4)i, 
(3.1.10),  (3.1.5),  (3.1.6)i  or  (3.1.6)2,  (3.1.11)  and  (3.1.9)i  or  (3.1.9)2.  Within  the 
context  of  this  investigation  it  will  be  assumed,  however,  that  V  or  is  chosen 
so  that  a  non-trivial  base  process  exists. 

3.2.  Perturbation  of  the  base  process.  Suppose  that  at  the  instant 
t  =  0  the  out-of-plane  displacement  and  velocity  fields  and  the  configuration  of  the 
phase  boimdaxy  associated  with  the  motion  specified  in  Section  3.1  are  subjected 
to  a  perturbation.  Let  this  perturbation  be  such  that  the  phase  boimdary  can 
be,  at  ^  =  O-f,  described  by  the  graph  Co  of  a  continuous  function  h  :  IR  IR 
of  the  X2“Coordinate,  and  segregates  elliptic  phases  of  the  three-phase  material 
at  hand  in  a  sense  consistent  with  that  which  was  present  for  t  in  (  — oo,0).  Let 
the  out-of-plane  displacement  and  velocity  fields  linked  to  this  perturbation  be 
given,  respectively,  by  a  once  continuously  differentiable  function  77  :  ]R^  IR 
and  a  continuous  fvmction  w  :  IR?  JR.  Assume  that  /i,  77  and  ^  represent 
small  deviations,  in  some  appropriate  sense,  from  their  counterparts  in  the  base 
process.  In  particular,  suppose  that  /i,  77,  77,0,  and  w  are  all  square  integrable 
on  their  domains  of  definition.  Furthermore,  require  that  the  components  of  the 
gradient  of  77  allow  the  satisfaction  of 

lim  77,0  (xi, 3:2 )^,a(a:i,  0:2)  =  0, 

Il+Xj— *00 

while  vj  complies  with 

lim  tn(xi,X2)  =  0,  (3.2.2) 

+  *00 

SO  that  the  disturbance  is  localized  in  a  neighborhood  of  the  phase  boundary 
associated  with  the  base  state  at  t  =  0. 

The  perturbation  at  <  =  0  will  initiate  a  new  process  involving  an  out-of- 
plane  displacement  field  u  :  IR^  x  1R+  — ♦  IR  which  is,  in  general,  a  weak  solution 
of  the  field  equation  (2.4.2)  £ind  satisfies  the  jump  conditions  in  (2.4.4)  and  (2.4.5) 
at  all  discontinuities  in  its  gradient,  the  kinetic  relation  (2.3.8)  or  (2.3.9)  at  sill 


(3.2.1) 
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pheise  boundaries,  and  the  initial  conditions 


ix(-, -,0+)  =  tio(-,0+)  +  7/  on 

•,  0+)  =  uo(",  0+)  +  on  R^ . 


(3.2.3) 


Since  the  perturbation  is  small,  assume  that  the  subsequent  process  involves 
only  a  single  phase  boundary  Ct  =  {(ri,r2,<)|  ^2  €  R}  for  eaeh  t  in 

R^ ,  with  <;  :  Rx  R+  — »  R  continuously  differentiable  on  its  domain  of  definition 
and  defined  so  that  it  complies  with  the  initial  condition 


<;(-,0+)  =  /i  on  R. 


(3.2.4) 


With  the  intent  of  linearizing  the  field  equation  in  (2.4.2)  about  the  base  process, 
write,  for  each  t  in  R-^-, 


u{xi,X2,t)  =  uo(xi,t)  +  w{xi,X2,t)  V(ii,X2)  €  P  \  Ct,  (3.2.5) 


where  w  zmd  its  derivatives  are  assumed  to  represent  small  departures  from  the 
relevent  quantities  in  the  base  process.  Assume  that  the  components  of  the  gra¬ 
dient  of  w  satisfy  the  following  limits: 

lim  lu,!  (xi, •, •)  =  0  on  RxR+, 

II  — *±00 

lim  w,2  (-,  X2,  •)  =  0  on  RxR+. 

12— ♦±00 

From  (3.2.3)  and  (3.2.5)  it  is  clear,  moreover,  that — when  inertial  effects  are  not 
ignored — the  increment  w  to  the  out-of-plane  displacement  field  must  satisfy  the 
following  initial  conditions: 


(3.2.6) 


w(-,-,0+)  =  7j  on  R^, 
=  tv  on  R^. 


(3.2.7) 
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It  is  important  to  emphasize  that  these  can  not  be  imposed  in  the  inertia-free 
setting. 

Next,  define  s  :  IRx  — >  iR,  the  correction  to  the  interface  position  due 

to  the  perturbation,  via 


=  vot  +  on  IR  Vt  €  iR+.  (3.2.8) 

Note,  from  (3.2.4)  that  the  increment  s  to  the  phase  boundary  position  must 
satisfy  the  initial  condition 


s{-,Q+)  =  h  on  IR.  (3.2.9) 

Observe  that  the  unit  normal  vectors  n±(-,t)  :  IR  M  to  Ct  are  given  by 
n±(-,()  =  ±^Tii==  on  It  V(€JR+.  (3.2.10) 

-f  5,2  (*,1) 

For  the  remainder  of  this  work,  choose  the  unit  normal  vector  associated  with  the 
plus  sign  in  (3.2.10)  emd  drop  this  sign  when  referring  to  it.  The  normal  velocity 
Vn{-,t)  :  ]R  —*  JR  of  Ct  is  given,  for  each  t  in  IR+,  by 

V-„(.,()  = on  M  V(€H+.  (3.2.11) 

3.3.  Linearization  of  the  field  equations  associated  with  the  process 
initiated  by  the  perturbation.  Let  V[  amd  PJ"  denote,  for  each  t  in  JR+,  plane 
sets  defined  by 


D'  =  {(i,,I2)|i,  <«(!,,<)),  V;  =  R^\Vl.  (3.3.1) 


Let  and  A'r  be  given,  in  turn,  by 


=  {(a:i,X2,<)|  (xi,X2)  €  V[,t  G  JR+}, 


(3.3.2) 
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and 

<^r  =  {(ij,X2,<)!  {^1,^2)  e  Vt,t  e  iR+}.  (3.3.3) 

o  o 

The  field  equations  which  hold  on  Xi  and  can  be  obtained  by  linearizing 
the  partial  differential  equation  in  (2.4.2)  about  7;  and  7r,  respectively.  First, 

O 

consider  the  derivation  of  the  field  equation  which  holds  on  Xi.  From  (2.4.3), 
(3.2.5)  and  the  assumption  regarding  the  magnitude  of  the  spatial  gradient  of  the 
increment  w  to  the  out-of-plane  displacement  field  it  is  clear  that 

7  =  -h  iw,i  )2  -I-  w,l  =  yj'yf  +  2yiw,i  +W,aW,a 

=  'Yi  +  w,i  on  Xi.  (3.3.4) 

From  (2.2.8),  (3.3.4)  and  Taylor’s  theorem  it  is  fiu-ther  evident  that 

M{^)  =  M{'yi  +  vj,i)  =  M{yi)  +  on  Xi.  (3.3.5) 

Next,  using  (3.2.5)  and  (3.3.5)  in  the  left-hand-side  of  the  partial  differential 
equation  in  (2.4.2)  gives 

),a  =  +  ■W'(7,)U>,1  )U,Q  ],o 

-  [(-^(tO  +  ){'ri  -I-  W,i  )],i 

-I-  [(M(7/)  +  )w,2  ],2 

=  r'(7j)w,ii  -|-M(7/)u;,22  on  Xi  (3.3.6) 

Note  that,  in  deriving  (3.3.6),  the  smoothness  of  r  and  hence  M,  the  identity 

r(7)  =  M(7)7  e  ]R+, 

which  follows  from  (2.2.6)  and  (2.2.8)  zind  its  consequence 


t'(j)  =  M(j)  +  M'{j)y  V7  G  \  {7,7} 
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have  been  used.  Observe,  also,  that  t'(7/)  >  0  by  whichever  of  (3.1.2)  is  appro¬ 
priate,  and  >  0  by  (2.2.9).  Prom  (3.1.1),  (3.2.5)  and  (3.3.6),  the  linearized 

O 

field  equation  which  holds  on  A’ j  is 


a]w,n+bjw,22=  w,  (3.3.7) 

where  the  positive  constants  a/  and  6/  are  defined  by 

ai  =  y/T'(-yi)/p,  bi  =  (3.3.8) 

O 

Similarly,  the  linearized  field  equation  which  holds  on  Ar  is 


ai.w,n  +K'^,22  =  w. 


(3.3.9) 


where  the  positive  constants  and  6r  are  defined  by 

a,  =  Vr'M/f,  6r  =  Vmyr)/p.  (3.3.10) 

In  writing  (3.3.10),  the  positivity  of  T'{yr)  >  0  and  M{yr)  >  0,  which  Eire  results 
of  (2.2.9)  and  whichever  of  (3.1.2)i,2  is  appropriate,  have  been  used. 

From  (2.4.7)  and  (3.3.6)  it  is  clear  that,  in  the  inertia-free  setting,  equations 
(3.3.7)  and  (3.3.9)  are  supplanted  by 


~i'bjw,22  —  0,  (3.3.11) 

and 

+blw,22  =  0,  (3.3.12) 

o  o 

which  hold  on  A/  and  Ar,  repectively. 

3.4.  Linearization  of  the  jump  conditions  and  kinetic  relation  as¬ 
sociated  with  the  process  initiated  by  the  perturbation.  Since  the  set 


-27- 


r  =  {(xi,X2,<)|  (a:i,a:2)  €  Ct,t  €  -K+l  represents  the  post-disturbance  trajec¬ 
tory  of  the  phase  boundary,  the  jump  conditions  in  (2.4.4)  and  (2.4.5)  and  the 
kinetic  balance  equation  in  (2.3.8)  or  (2.3.9) — with  Vn  and  /  given,  respectively, 
by  (3.2.11)  and  (2.4.6) — must  hold  on  it.  Assume,  henceforth,  that  the  function 
s  introduced  via  (3.2.8)  and  its  derivatives  are  small  in  the  same  sense  that  w  is 
small.  Note,  first,  that  this  assumption  implies,  using  (3.2.10)  and  (3.2.11),  the 
following  approximations  for  n  and  V„  on  F : 

n  =  61—5,2  62  on  r,  Vn  =  +  •s  on  r.  (3.4.1) 

It  will  now  be  shown  that  a  further  consequence  of  the  above  stipulation 
regarding  the  size  of  s  and  its  derivatives  is  that,  within  the  error  associated  with 
the  linearization,  the  jump  conditions  in  (2.4.4)  and  the  kinetic  relation  can  be 
enforced  on  an  undisturbed  continuation  of  the  phase  boundary  intrinsic  to  the 
base  process  I  given  by 

/  =  {(xi,X2,t)|xi  =  Vot,X2  €  G  JR+}-,  (3.4.2) 

but,  when  Vn  appears  in  any  of  these,  the  contribution  due  to  s  from  (3.4. 1)2 
must  be  retained.  To  see  this  consider,  for  example,  the  limiting  values  of  the 
xi-component  of  the  gradient  of  the  out-of-plane  displacement  field  u  on  either 
side  of  the  phase  boundary;  note  that  by  Taylor’s  theorem,  (3.1.1),  and  (3.2.8), 

u,i  (<;(x2,  t)-,  X2,  i)  :=  lim  u,i  (t>o<  -  h  +  5(12,  t),  X2  -t-  hs,2  (x2,  t),  t) 

h\,0 

=  7/  +  ((^oi  -1- 

=  71  +  (vot-,X2,t)  -t-  rv,ll  (voi-,X2,t)s(x2,0 

=  U,1  (vot-,X2,i)  +  W,11  {vot-,X2,t)s{x2,t) 

V(x2,t)  € -K  X  (3.4.3) 


-28- 


and,  similaxly, 

«,i  X2,t)  :=  lim  u,i  {vot  +  h  +  s{x2,t),  X2  +  hs,2  (x2,t),  t) 

h\,0 

=  7r+  ((vot  +  -5(^2,  ^2,0 

=  7r  +  W,i  ivot+,X2,^)  +  If, 11  ivot  +  ,X2.,t)s{x2,t) 

=  11,1  {vot+,X2,t)  +  U>,n  {vot  +  ,X2,t)s{x2,t) 

V{x2,i)e  Rx  M+.  (3.4.4) 

Since  both  s  and  the  derivatives  of  w  are  assumed  smzJl  and  of  the  same  order,  the 
quadratic  terms  on  the  right-hand-sides  of  (3.4.3)  and  (3.4.4)  can  be  neglected. 
This  produces 

u,i{<;{x2,t)±,X2,t)  =  u,i  {vot±,X2,t)  W{x2,t)  6  R  X  R^.  (3.4.5) 

Within  the  scope  of  the  linearization,  u,i  {<;{x2,t)dt,X2,t)  is,  therefore,  obtained, 
for  each  (x2,t)  in  Rx  JR+,  by  evaluating  u,i  (•,X2,t)  at  uot±.  Analogous  remarks 
also  hold  for  the  r2“Component  of  the  gradient  of  the  out-of-pleine  displacement 
field,  the  out-of-plane  velocity  field  and  the  shear  strain  field.  Hence,  with  the 
aid  of  (3.4.1),  (2.4.4)  implies  that 

0  =r(7r)  -  r(7i)  -f  pvo(vr  -  vi)  +  pvo{jr  -  7i)s{x2,t) 

+  p{alw,i  (vot+,X2,t)  -  ajw,i  {vot-,X2,t)) 

+  pvo{w(vot+jX2,t)  -  w(vQi-,X2,t))  V{x2,t)  e  R  X  R^.  (3.4.6) 

From  (3. 1.4)2,  the  constant  term  on  the  right-hand-side  of  (3.4.6)  is  zero  and, 
hence,  the  linearization  of  the  jump  condition  which  enforces  the  balance  of  Unear 
momentum  across  the  phase  boundary  leads  to 

0  =alw,i  {vQt+,-,t)  -  afw,i  (vo<-,i2,0  +  vo(7r  -  7l)s{^2,t) 


+  Vo{w{vot  +  ,X2,t)  — 'U’{vot  —  ,X2,t))  V(x2 ,  t)  G  iR  X  .  (3.4.7) 
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Prior  to  deriving  the  linearized  kinetic  relation  it  is  convenient  linearize  the 
driving  traction  /.  From  (3.2.4),  (2.4.3)  and  the  foregoing  discussion  one  finds 
that 

f{x2,t)  =  fo  +  \p{v  -  lT){alw,i  {vot-\-,X2,t)  +  ajw,i  {vot-,X2,t)) 

+  2(^(>)  ~  ’■(7/))(«>U  {vot+,X2,t)  -  W,1  (vot-,X2,t)) 

V(x2,t)  e  Mx  JR+,  (3.4.8) 

where  the  base  driving  traction  /o  is  given  by  (3.1.7).  From  (3. 1.4)2  it  is  clear, 
furthermore,  that  (3.4.8)  simplifies  to  read 

f(x2,t)  =  fo  +  ^P(7I  -  7r)((ar  “  (^0*+, 3:2,0  +  (of  -  Vo)w,i  (vot-,X2,t)) 

V(x2,t)  6  JRx  JR+.  (3.4.9) 

If  the  kinetic  relation  is  of  the  form  given  in  (2.3.8),  then  (3.4. 1)2,  (3.4.9) 
and  Taylor’s  theorem  lead  to 

Vo  +  5(3^2, 0  -y(fo) 

+  \pV'{fo){7i  -  lr)({al  -  vl)w,i  {vot+,X2,t)  +  (af  -  vl)w,i  {vot-,X2,t)) 

y{x2,t)e  Mx  1R+.  (3.4.10) 

If,  on  the  other  hand,  the  kinetic  relation  is  provided  in  the  form  (2.3.9),  then 
(3. 4. 1)2,  (3.4.9)  and  Taylor’s  theorem  give,  simileirly, 

<p(vo)  +  ‘f'(vo)s{x2.,t)  =  /o 

+  ^p(7l  -  7r){(al  -  UoHi  (vot+, 3:2,0  +  (of  -  u2)«’n  (uot-, 3:2,0) 


y(x2,t)  e  Six  Si+. 


(3.4.11) 
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Use  of  the  (3. 1.6)1  and  (3. 1.6)2  in  (3.4.10)  and  (3.4.11),  respectively,  results  in 
the  linearized  kinetic  relation 


i(x2,0 


7<  -  7r 

2t;* 


{{al  -  Vo)w,i  (vot+,X2,t)  +  (o/  -  Vo)w,i  (vot-,X2,t)) 


V(x2,t)  G  M  X  JR+, 


(3.4.12) 


where  the  constant  v»  is  defined  by  either 

_  1 

''*"pU'(/o)’ 

if  the  kinetic  relation  is  furnished  in  the  form  (2.3.8),  or 


(3.4.13) 


P 


(3.4.14) 


if  the  kinetic  relation  is  supplied  in  the  form  (2.3.9).  Note,  from  (3.1.9),  that  v, 
is  a  real — but  nonzero — constant. 

Consider,  now,  the  task  of  linearizing  the  jump  condition  in  (2.4.5).  Note, 
from  (3.2.5),  that  (2.4.5)  implies 


0  =uo(C(i2,0+,t)  -  uo(((x2,t)-',i)  +  «^(C(a:2,0+> 3:2,0  - 

V(x2,t)e  JRx  JR+,  (3.4.15) 

Certainly,  from  (3.1.1),  (3.2.8)  and  (3.1.4)i, 

«o(C(3;2,0+,0-“o(C(3^2,0-,0  =  (7J-7r)5(x2,0  V(i2,0  E  JRxJR+;  (3.4.16) 

furthermore,  from  (3.2.8),  Taylor’s  theorem,  and  the  assumption  regarding  the 
small  magnitude  of  products  involving  s  and  the  derivatives  of  w, 

w((;(x2,t)±,X2,t)  =  w(vot±,X2,t)  +  (vot±,X2,t)s(x2,t) 


=  w(vot±,X2,t)  V(x2,t)  €  JR  X  JR^. 


(3.4.17) 
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Hence,  from  (3.4.15)-(3.4.17),  the  linearization  of  the  jump  condition  which  en¬ 
forces  the  continuity  of  displacement  across  the  phase  boundary  yields 

'w{vot  +  ,X2,t)  - ‘w{vot-,X2,t)  =  ('fl  - ‘yr)six2,t)  ^{x2,t)  E  ]R  X  1R+.  (3.4.18) 

Differentiation  of  (3.4.18)  with  respect  to  time  then  results  in  the  following 
identity: 


0  =w(vQt+,X2,t)  -  w{vot-.,X2,t)  +  (7r  “  'yi)s(x2,t) 

+  Vq{w,i  {vQt  +  ,X2,t)  -  W,1  (vot-,X2,t)) 

V(x2,t)  e  JR  X  JR+,  (3.4.19) 

Appropriate  substitution  of  (3.4.19)  into  the  linearization  of  the  jvunp  condition 
which  enforces  the  balance  of  linear  momentum  across  the  phase  boundary  (3.4.7) 
then  gives  rise  to 


0  =(a^  -  v^)w,i  (vot+,X2,t)  -  (af  -  t;o)u;,i  (vot-,X2,t) 

-t- 2uo(7(  -  7r)s(i2,0  V(r2,t)  €  iR  X  iR+.  (3.4.20) 

By  virtue  of  the  foregoing  calculations  it  is  crucial  to  note  that,  within  the 
scope  of  the  linearization,  it  is  legitimate  to  enforce  the  partial  differential  equa¬ 
tions  in  (3.3.7)  and  (3.3.9)  on  the  interiors  of  the  sets  1?/  eind  i?r  defined  by 

12/  =  {(xi,Z2,<)|(xi,r2)  E  Jli,i  e  JR+},  (3.4.21) 

with  JTj  =  {(ii,Z2)lxi  <  Vot,X2  €  JR)  for  each  t  in  JR^.,  and 

12r  =  {(zi,Z2,t)l(zi,Z2)  6  ni,t  €  iR+},  (3.4.22) 

O 

with  ni  =  {(zi,Z2)lzi  >  Vot,X2  €  iR}  for  each  t  in  iR+,  instead  of  the  sets  Xi 
and  Xr  (recall  (3.3.2)  and  (3.3.3)). 
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In  the  inertia-free  case  it  is  readily  shown  that,  while  (3.4.18)  continues  to 
hold,  (3.4.17)  is  replaced  by 


a\w^i{vQi+^X2,t)  —  aiW,i{vQt—,X2,t)  =  Q  V(x2,t)  E  M  x  ]R+,  (3.4.23) 


and  (3.4.11)  simplifies  to  read 

s(x2,t)  (vot-^,X2.,t)  -I-  ajw,i  (vot-,X2,t)) 

y{x2j)e  IRx  1R+.  (3.4.24) 


Finally,  remarks  analogous  to  those  made  regarding  the  enforcement  of  the  partial 
differential  equations  in  (3.3.7)  and  (3.3.9)  on  1^/  and  hr  apply  also  to  those  in 
(3.3.11)  and  (3.3.12). 

3.5.  Linearized  description  of  the  post  perturbation  process.  In 
this  section  the  linearized  field  equations,  jump  conditions,  kinetic  relation,  ini¬ 
tial  conditions  (where  appropriate),  and  far  field  decay  conditions  satisfied  by 
the  increments  w  and  5  to  the  out-of-plane  displacement  field  and  the  interface 
position  eire  listed  in  both  the  inertizJ  and  inertia-free  cases. 

In  the  inertial  case,  (3.3.7)  and  (3.3.9)  give  the  following  linearized  field 
equations 

+1»/W,22  =  ^  on  1^/, 

(3.5.1) 

al'w,u  +blw,22  =  w  on  l?r. 

In  addition,  from  (3.4.20)  and  (3.4.18),  the  following  jump  conditions  hold 


-  Uo)u^,i  I  =  2t;o(7r  -  7/)^  on  /, 
[wj  =  (7,  -  7r)s  on  /, 


(3.5.2) 


where 


a 


2 

;  • 


(3.5.3) 
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Next,  from  (3.4.12)  the  following  linearized  kinetic  relation  holds: 

s  =  — — ))  on  I.  (3.5.4) 

The  initial  conditions  satisfied  by  w  and  s  are,  from  (3.2.4)  eind  (3.2.7), 

io(-, -,0+)  =  t;  on  IR^, 

iw(-, -,0+)  =  07  on  (3.5.5) 

s(-,0+)  =  h  on  M. 


Finally,  from  (3.2.6),  it  is  assumed  that  the  following  far  field  conditions  hold 


hm  10,1  (xi,  •,<)  =  0  on  iR, 

II  —>±00 

lim  to,2  (•,X2,<)  =  0  on  iR, 

ij— »±oo 

for  each  t  in  JR+. 

In  the  inertia-free  case,  (3.5.1)  is  replaced  by  the  following: 

<ifu>,ii+bfw,22=0  on 
0^10,11  +6^10,22  =  0  on  j^r. 


(3.5.6) 


(3.5.7) 


Furthermore,  the  jump  condition  (3.5.2)]  is,  by  virtue  of  (3.4.23),  replaced  by 


[a^io,i  1=0  on  /, 


(3.5.8) 


while  (3. 5. 2)2  continues  to  hold.  Following  (3.4.24),  the  linearized  kinetic  relation 
(3.5.4)  is  superseded  by 


5  =  21 — on  /,  (3.5.9) 

In  the  absence  of  inertial  effects  initial  conditions  cannot  be  given  for  the  incre¬ 
ments  to  the  out-of-plane  displacement  and  velocity  fields  w  and  w\  the  initial 


condition  (3. 5. 6)3  pertaining  to  s  still,  however,  continues  to  be  applicable.  The 
decay  conditions  (3.5.6)  also  still  hold. 

3.6.  Normal  mode  analysis  in  the  inertia-free  setting.  An  approx¬ 
imate  means  for  analyzing  the  linear  stability  of  the  base  process  described  in 
Section  3.1  is  afforded  by  the  study  of  the  inertia-free  initial  value  problem  con¬ 
sisting  of  (3.5.7)-(3.5.9),  (3. 5.2)2,  (3.5.5)3  and  (3.5.6).  Observe  that,  by  virtue 
of  the  linearization,  the  relevant  partial  differential  equations,  jump  conditions 

and  kinetic  relation  are  cill  linear  with  constant  coefficients;  note,  also,  that  the 
0.0 

domains  J7J  and  are,  for  each  t  in  M+,  rectangular.  It  is  therefore  possible 
to  find  a  solution  to  the  linearized  pjirtial  differential  equations,  jump  conditions 
and  kinetic  relation  in  the  form 


,  ,  f  V(ii,X2)€^i,  t  G  M+, 

w(Xl,X2,t)  =  <  ,  ,  ,  •  o 

V(xi,X2)€/7r,  t  6  M+,  (3.6.1) 

s(x2,  <)  =  \/{x2,t)  €  IR  X  JR+, 

where  the  amplitudes  Wi,  Wr  and  S',  wave-numbers  and  k,  and  growth-rate 

p  are  all  constants.  To  comply  with  the  decay  condition  (3. 5. 6)1  it  is  clear  that 
3fJ(^f)  and  3?(^r)  must  both  be  positive.  The  Ansatz  (3.6.1)  is  not,  in  general, 
consistent  with  the  initial  condition  (3.5. 5)3  or  the  decay  condition  (3. 5. 6)2;  since 
the  initial  disturbance  h  is  stipulated  to  be  square  integrable  on  iR,  and  hence 
can  be  represented  as  a  Fourier  integral — 


+  00 

h{x2)  =  TT  /  dK  Vx2  e  iR, 

2n  J 


(3.6.2) 


it  is  reasonable  to  expect  that  stability  results  cem  be  obtained  by  a  normal¬ 
mode  analysis;  such  an  analysis  entails  substitution  of  (3.6.1)  into  (3.5.7)-(3.5.9) 
to  determine  the  growth-rate  p  as  a  fimction  of  the  positive  wave-numbers 
and  and  the  real  wave-number  k.  In  the  context  of  such  an  undertaking,  the 
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amplitude  S  and  wave-number  k  axe  regarded  as  given  and  non-zero,  while,  along 
with  the  growth-rate  p,  the  amplitudes  Wi  and  Wr  and  wave-numbers  and 
are — due  to  the  present  lack  of  inerticil  effects — to  be  determined.  If  there  exists 
a  complex  growth-rate  p  with  positive  real  part  which  arises  els  a  solution  to 
the  aforementioned  problem  then  the  base  process  will  be  referred  to  els  linearly 
unstable.  Otherwise,  the  base  process  will  be  called  linearly  stable. 

Substitution  of  (3.6.1)  into  (3.5.7)-(3.5.9)  and  (3. 5. 2)2  yields  the  following 
system  of  five  equations  in  the  five  unknowns  IV; ,  Wr,  ir  and  p: 

(aUf  -  bfK‘^)W,  =  0, 

(aUJ  -  bW)Wr  =  0, 

«rCrfVr  -|-  a](iWl  =  0,  (3.6.3) 

Wr-W,-  (7,  -  7;)5  =  0, 

-  aUlWi)  +  Sp  =  0. 

First,  (3.6.3)i,2  give,  recalling  that  and  must  be  positive, 

6  =  — l«l,  =  — |k|.  (3.6.4) 

Ql 


In  particular,  (3.6.4)  implies  that,  the  wave- numbers  and  must  be  real  and, 
further,  that,  of  and  k,  the  growth-rate  p  will  depend  only  on  k.  Next, 

(3. 6. 3)3, 4, 5  can  be  solved  for  the  remaining  unknowns  W;,  IVr  and  p  to  yield 


where  the  constant  is  defined  as  follows: 


^2  _  -  Irf 

aibi  -|-  ttrbr 


(3.6.6) 
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Equation  (3. 6. 5)3  which  gives  the  growth-rate  p  in  terms  of  the  wave-number 
K  and  various  physical  parameters  intrinsic  to  the  problem  at  hand  will  be  referred 
to  as  the  dispersion  relation.  It  is  clear  from  (3. 6. 6)3  that  the  growth-rate  p  is 
real.  Since  k  ^  0  by  assumption  and  /  0  by  (3.1.10),  it  is  clear,  moreover, 
that  the  signs  of  p  and  u,  are  opposite:  viz., 

sgn(p)  =  sgn(-u,).  (3.6.7) 

The  linear  stability  of  the  base  process,  in  the  absence  of  inertial  effects,  thus 
depends  entirely  upon  the  sign  of  u,.  Significantly,  the  wave-number  k  plays  no 
role  in  determining  stability.  Moreover,  the  local  mechanical  properties  of  the  high 
and  low  strain  elliptic  phases  of  the  three-phase  material  at  hand  do  not  effect  the 
stability  criteria.  If  u,  >  0  then  the  base  process  is  linearly  stable  with  respect  to 
all  initial  disturbances  h  of  the  type  under  consideration.  If,  alternatively,  u*  <  0 
then  the  base  process  is  linearly  unstable  with  respect  to  all  disturbances  h  of 
the  type  considered  here.  To  summeirize,  in  the  absence  of  inertial  effects,  the 
criterion  u,  <  0  is  necessary  and  sufficient  for  the  base  process  to  be  unstable 
with  respect  to  any  perturbation  of  the  type  put  into  consideration  in  Section 
3.2.  Thus,  if  the  kinetic  response  fimction  V  or  ip — as  is  appropriate  to  whether 
a  kinetic  relation  of  the  form  (2.3.8)  or  (2.3.9)  is  prescribed — is  non-monotonic, 
the  base  process  may  be  linearly  unstable  with  respect  to  any  initial  disturbance 
h  of  the  type  being  considered.  Discussion  regarding  the  physical  reasonableness 
of  a  non-monotonic  kinetic  response  fimction  is  left  until  Section  3.9. 

Suppose,  now,  that  the  normal  velocity  uq  of  the  phase  boundary  in  the  base 
process  is  zero.  Note  that  this  is  equivalent  to  requiring  that  the  base  process 
be  mechanically  equilibrated.  Then,  by  the  admissibility  of  the  kinetic  response 
function,  one  or  both  of  fo  and  u,  must  be  zero.  Since  the  latter  contradicts 
(3.1.10),  uo  =  0  implies,  at  present,  that  /o  =  0  and,  furthermore,  either  ^”'(0)  >  0 
or  >  0.  Thus  the  foregoing  stability  dichotomy  implies  that  a  mechanically 

equilibrated  base  process  of  the  type  defined  in  Section  3.1  must  be  linearly  stable 


-37- 


with  respect  to  all  initial  disturbances  h  of  the  type  considered  heie. 

Note  that  the  foregoing  results  are  consistent  with  those  presented  by  Fried 
[14]  in  a  study  of  the  distribution  of  driving  traction  along  a  particular  non- 
planax  phase  boimdary.  The  latter  work  shows,  roughly  speaking,  that  the  driving 
traction  is  less  in  regions  of  the  interface  where  the  curvature  is  larger.  Specifically, 
if  small  deviations  of  the  relevant  type  from  a  planar  interface  with  constant  base 
driving  traction  /*  are  considered  and  a  kinetic  response  function  V  is  provided 
so  that  <  0  then  the  normal  velocity  of  the  interface  in  regions  of  higher 

curvatme  exceeds  that  in  regions  of  lesser  curvature.  The  interface,  then,  has  a 
marked  tendency  to  evolve  in  a  manner  wherein  its  curved  portions  move  ahead 
of  its  flat  portions — and,  hence,  become  less  planar.  Such  a  response  is  intuitively 
unstable.  K,  however,  /,  is  such  that  y'(/*)  >  0,  flat  portions  of  the  interface 
tend  to  catch  up  with  the  ciurved  portions  so  that  the  interface  regains  its  planar 
shape — hence,  the  planar  interface  is  stable. 

3.7.  An  alternative  to  normal  mode  analysis  in  the  inertia-free 
setting.  The  analysis  performed  in  Section  3.6  resulted  in  a  necesary  and  suffi¬ 
cient  condition  for  the  base  process  to  be  unstable  with  respect  to  an  arbitrary 
perturbation  of  the  type  discussed  in  Section  3.2;  it  did  not,  however,  encom¬ 
pass  a  means  for  tracking  the  evolution  of  either  a  stable  or  imstable  response  to 
perturbation.  This  section  will  focus  on  an  analysis  which  does  allow  the  post¬ 
perturbation  evolution  of  the  phase  boundary  to  be  followed  in  the  linear  regime. 
Since  the  analysis  is  performed  in  the  linear  realm  it  is,  in  the  case  of  an  imsta¬ 
ble  response,  only  of  value  in  a  short  time  intervEil  following  the  perturbation  at 
t  =  0.  Consider,  agadn,  the  inertia-free  initial  boundary  value  problem  consisting 
of  (3.5.7)-(3.5.9),  (3.5. 2)2,  (3. 5. 5)3  and  (3.5.6).  Techniques  from  potential  the¬ 
ory  arc  used  in  the  appendix  of  this  work  to  show  that  the  increment  w{-,-,t) 
to  the  out-of-plane  displacement  field  can  be  represented,  for  each  t  in  iR+,  in 
terms  of  the  sum  of  a  single-layer  potential,  with  density  proportionzil  to  i,  and  a 
double-layer  potential  potential,  with  density  proportional  to  s — each  on  the  line 
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x\  =  vot — as  follows: 


+  00 

w{xi,X2,t)=  J  Ki(xi  -VQt,X2  - 


t)dC 


—  oo 
-foo 


+  J  K2(xi  -vot,X2  -  '^{xi,X2)  e  nf.{Z.7.1) 

—  oo 

The  sub-  and  superscripts  d  in  (3.7.1)  are  to  be  replaced  by  either  /  or  r — as 
appropriate,  while  the  kernels  Ka{-y’,c)  :  IR^  \  {(0,0)}  — +  M  are  defined  via 

^  Inv/c2x2-hy2  V(x,y)eiR2\{(0,0)}, 

(3.7.2) 


A'i(x,  y,c)  = 


K2{x,y,c)  = 


2ir  1/2(7,  -  7r) 

1  cx 


27r  c2x2  ^  y2 

for  each  c  in  ]R+,  and  s,  introduced  in  (3.2.8),  satisfies  the  following  functional 
initial  value  problem 


(3.7.3) 


2 

i(x2,0  =  —  \/{x2,t)eJRxR+, 

TTU*  J  X2  —  C 

—  oo 

s(l2,  0)  =  /l(x2)  VX2  e  iR. 

The  integral  on  the  right-hand-side  of  (3.7.3)i  is,  as  indicated,  of  the  Cauchy 
principal  value  type. 

The  stability  of  the  base  process  rests  on  the  stability  of  the  solution  s  to 
the  foregoing  functional  initial  value  problem.  To  investigate  this  issue  assume — 
because  of  the  square  integrability  of  h — that  s  can  be  represented  in  the  form 

+  O0 

=  f  s(K,<)e’^dK  V(x2,t)  6  ^  X  .K+,  (3.7.4) 

27r  J 


for  some  fimction  s  :  IR  x  JR+  — »  M  which  satisfies 

+  00 

s(/c,0)  =  J  /i(x2)e~''‘**  dx2  =:  h{K)  V/c  €  IR. 


(3.7.5) 
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Now,  from  (3.7.3)-(3.7.5)  it  is  possible  to  derive,  for  each  fixed  k  in  JR,  the 
following  initicd  value  problem: 


i(/c,-)H - jK|5(/c,-)  =  0  on  JR, 

u* 

s(k,0)  =  h{K). 


(3.7.6) 


Observe  that  the  Laplace  transform  could  be  applied  to  (3.7.6)  to  derive  the  dis¬ 
persion  relation  (3.6.5)3.  The  linear  stability  results  obtained  in  Section  3.6  would, 
then,  follow  immediately.  Instead,  note  that  (3.7.6)  can  be  solved,  formally,  to 
yield  the  following  expression  for  s: 


J(/c,<)  =  h(/c)e-^l'‘l‘  y{K,t)e  IRx  1R+. 


(3.7.7) 


Inspection  of  (3.7.7)  reveals  that,  as  demonstrated  in  Section  3.6,  the  linear  sta¬ 
bility  of  the  base  process  is  decided  entirely  by  the  sign  of  u,.  Note  that  if  y,  >  0 
then,  by  (3.7.4)  and  (3.7.7), 


K^2,<)  =  —  / 

Try,  J 


HQdC 

i^2-cy+i^y 


V(x2,0  €  .K  X  JR+. 


(3.7.8) 


Hence,  under  these  circumstances,  it  is  clear  that 


lim  s{x2,t)  =  0  Vi2  G  JR, 

t— ‘OO 


(3.7.9) 


for  all  square  integrable  initial  disturbances  h.  The  base  process  is,  therefore, 
linearly  stable  with  respect  to  all  such  disturbances. 

Suppose,  now,  that  y,  <  0 — which,  as  remarked  in  Section  3.6,  can  occur 
only  if  Vo  >  0.  It  is  then  instructive  to  consider  a  special  example  where  the 
initial  disturbzmce  h  is  given  by 


(3.7.10) 
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with  t  taken,  without  loss  of  generality,  to  be  a  positive  constant.  See  Figure  2 
for  a  representative  graph  of  h.  Note  that  since  (3.7.3)i  is  linear  it  is,  in  general, 
possible  to  consider  an  initial  disturbance  which  involves  a  linear  superposition 
of  disturbances  of  the  type  defined  in  (3.7.10),  e.g.,  H  •.  IR  —*  IR  given  by 


Vr2  ^  JR, 


where  iV  is  a  natural  number  and  A„,  y„  and  are  real  constants  for  each  n 
in  {1,2,...,  N}.  The  results  which  follow  are  easily  extended  to  apply  to  such 
a  generalization  of  the  initial  disturbance.  Now,  the  Fourier  transform  h  of  the 
function  h  defined  in  (3.7.10)  is  given  by 


hiK)  =  V/ceiR. 


(3.7.11) 


Next,  (3.7.4),  (3.7.7)  and  (3.7.11)  imply  that 


+  O0 

=  ^  V(X2,<)  e  iRx  [0,<c),  (3.7.12) 

27r  J 


where  the  critical  time  tc  is  chosen  small  enough  so  that  the  amplitude  of  s  does 
not  contradict  the  assumptions  necessary  for  the  linearization  to  remain  valid  over 
the  time  interval  [0,tc);  *-e., 


|v«|£ 


(3.7.13) 


From  (3.7.12)  it  transpires,  further,  that 


(3.7.14) 


The  expression  for  on  JR  given  in  (3.7.14)  shows  that  the  amplitude 

and  wavelength  of  the  hump  associated  with  the  special  initial  disturbance  imder 
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consideration  increase  and  decrease,  respectively,  as  t  approaches  tc.  Hence,  up 
until  some  critical  time,  the  assumption  u,  <  0  implies  that  the  hump  associated 
with  the  initial  disturbance  grows  in  an  imstable  fashion.  See  Figure  4  for  a 
depiction  of  this  unstable  evolution.  Observe,  moreover,  that,  as  t  approaches 
too,  the  function  s  given  by  (3.7.14)  forms  a  delta  sequence  so  that 

lim  s(-,t)  =  on  IR,  (3.7.15) 


where  S  is  the  Dirac  delta  distribution. 

3.8.  Normal  mode  and  energy  analysis  with  inertial  effects  in¬ 
cluded.  The  applicability  of  the  stability  criterion  obtained  in  Section  3.6,  and 
recovered  in  Section  3.7,  is,  so  far,  limited  to  contexts  where  inertial  effects  are  in¬ 
significant.  This  section  is  concerned  with  the  extension,  as  far  as  possible,  of  the 
inertia-free  stability  dichotomy  based  upon  the  sign  of  u,  into  the  inertial  realm. 
As  a  first  step  toward  achieving  this  goal  a  normal  mode  analysis  analogous  to 
that  performed  in  Section  3.6  is  undertaken  with  inertial  terms  present.  The 
conclusion  of  this  analysis  is  that  v*  <  0  is  necessary  and  sufficient  to  guarantee 
the  linear  instability  of  the  base  process  with  respect  to  a  particular  subset  of  the 
class  of  perturbations  put  into  consideration  in  Section  3.2.  It  is  then  eurgued  that 
the  remaining  class  of  perturbations  not  covered  by  the  normal  mode  analysis  is, 
in  fact,  very  small.  Next  an  energy  argument  is  used  to  show  that  u,  <  0  is  a 
necessary  condition  for  the  base  process  to  be  linearly  unstable  with  respect  to  an 
arbitrary  disturbance  within  the  full  cleiss  of  perturbations  introduced  in  Section 
3.2. 

Consider,  now,  the  initial  boundary  value  problem  composed  by  (3.5.1), 
(3.5.2)  and  (3.5.4)-(3.5.6).  It  is,  once  again,  possible  to  find  a  solution  to  the 
system  formed  by  the  linearized  field  equations,  jump  conditions  and  kinetic  re¬ 
lation  in  the  form  of  (3.6.1).  Remarks  regarding  the  decomposition  of  the  initial 
data — which  now  includes  r]  and  za  as  well  as  h — and  the  satisfaction  of  (3. 5. 6)2 
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eikin  to  those  made  in  Section  3.6  are  pertinent.  Observe  that,  in  the  inertial  case, 
the  amplitudes  Wr,  Wi  and  S  and  wave-numbers  ^/,  and  k  must  be  viewed  as 
given  for  the  normal  mode  analysis  to  be  effective  in  determining  necessary  and 
sufficient  conditions — via  a  dispersion  relation  analogous  *o  (3. 6. 6)3 — for  insta¬ 
bility  with  respect  to  arbitrary  disturbances  of  the  out-of-plane  displacement  and 
velocity  fields  and  the  interface  within  the  class  of  perturbations  put  forth  in  Sec¬ 
tion  3.2.  Substitution,  however,  of  (3.6.1)  into  (3.5.1),  (3.5.2)  and  (3.5.4)  produces 
the  following  relations 


^  //(k,p)-uoP 

s/  —  2  2  ' 

af  - 

^  /r(K,p)  +  Vop 


w  =  (7<  -7r)(t>0p-/r(K,p))^ 

'  /K«,P)  +  /r(«,P)  ’ 

w  ^  (7<  -7r)(t;op-f /<(k,p)) 

//(«,p)  +  /r(«,p) 

(7/  -  7r)^(//(«,p)/r(«,p)  +  VqP'^) 

+  fr{K,p)) 


(3.8.1) 


where  and  fr  :  M  x  (E  —y  (T  eire  given  by 


/r(«,P) 


y/{a]  -v^)bfK^  +  a] 
yj {o-l  -  vl)hlK^ alp'^ 


V(«,p)  e  Mx  (E, 
V(k,p)  £  JRx  (E. 


(3.8.2) 


The  square  roots  in  (3.8.2)  are  defined  so  that  when  their  respective  arguments 
are  real  and  positive  the  resulting  square  root  is  also  real  and  positive. 

It  is  clear  from  (3.8.1)i,2,3,4  that,  for  (3.6.1)  to  represent  a  solution  to  (3.5.1), 
(3.5.2)  and  (3.5.4),  the  amplitudes  Wi  and  Wr  and  wave-numbers  and  cannot 
be  chosen  independently  of  5  and  k.  Hence,  the  normal  mode  procedure  is  only 
capable  of  analyzing  the  linear  stability  of  the  base  process  with  respect  to  a 
certain  class  of  perturbations;  that  is,  it  is  only  possible — via  this  analysis — to 
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determine  conditions  necessary  and  sufficent  for  the  instability  of  the  base  process 
with  respect  to  a  subset  of  the  class  of  perturbations  put  into  consideration  in 
Section  3.2.  To  do  this  it  suffices  to  einalyze  the  zero  structure  of  the  dispersion 
relation  (3. 8. 1)5  as  a  function  of  the  growth-rate  p  for  fixed  values  of  the  wave- 
niimber  k  and  the  pareuneters  7/,  7r,  uoi  o.l,  Od  K,  p  and  u*.  This  is  done 
below. 

K  u*  <  0  it  is  evident,  by  inspection,  that  there  exists  a  real  positive  root  p 
to  (3.8. 1)5  for  all  admissible  values  of  k,  7/,  7r,  uo,  a/,  Or,  bi,  br  and  p.  Hence, 
u,  <  0  is  sufficient — indepedent  of  the  value  of  the  wave-number  k — to  guarantee 
that  the  base  process  is  unstable  with  respect  to  the  n arrowed  class  of  initial 
disturbances  at  hand.  To  establish  the  converse,  show  that  the  condition  u*  >  0 
is  sufficent  to  guarantee  that  there  cannot  exist  unstable  zeros  p  to  (3.8.1  )5.  Let 
F(k,  •)  :  (T  — ♦  (T  be  given,  for  each  k  in  iR,  by 


F(k,p)  = 


(7<  -  7r)^(/<(p,  «)/r(p,  k)  -t-  ngp^) 
flip,  «)  +  /r(p,  «) 


Vp€  (T. 


(3.8.3) 


Now,  it  is  easy  to  show  that  if  3?(p)  >  0  then  3t(F(/c,p))  >  0  for  every  k  in  Ft. 
Hence,  if  u,  >  0  then  all  roots  p  to  (3.8. 1)5  must  have  non-positive  real  parts. 
Thus,  V,  <  0  is  also  a  necessary  condition  for  the  base  process  to  be  unstable 
with  respect  to  an  arbitrary  element  of  the  class  of  perturbations  which  can  be 
tested  via  the  normal  mode  analysis.  Observe  that  this  conclusion  is  independent 
of  the  wave- number  «  aissociated  with  (3.6.1). 

If,  in  place  of  the  foregoing  normal  mode  analysis,  a  full-fledged  Fourier- 
Laplace  transform  analysis  of  (3.5.1),  (3.5.2)  and  (3.5.4)-(3.5.6)  is  performed, 
then  the  narrowing  of  the  class  of  initial  data  necessitated  by  the  normad  mode 
analysis  does  not  occur.  Furthermore,  in  this  case  it  transpires  that  the  Fourier- 
Laplace  transform  of  s  can  be  expressed  in  the  form 


(3.8.4) 
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where  h  is  the  Fourier  transform  of  h  and,  for  each  (k,p)  in  iR  x  (T,  H{k,p)  is  a 
functional  of  the  initial  data  rj  and  w.  Evidently,  it  is  possible  that  there  exist 
combinations  of  p,  re?  and  h  which  would  allow  the  cancellation  of  an  unstable 
zero  in  the  denominator  of  the  expression  on  the  right-hand-side  of  (3.8.4)  by  a 
zero  in  its  numerator.  It  is  equally  clear,  however,  that  the  set  of  such  initial 
data  constitutes  a  very  small  one  within  the  full  class  of  initial  data  purveyed 
in  Section  3.2.  Hence,  except  in  response  to  very  special  initial  perturbations 
the  condition  u,  <  0  is  necessary  and  sufficent  for  the  base  state  to  be  linearly 
unstable.  The  normal  mode  analysis  thus  shows  that  u,  <  0  is  necessary  and 
sufficient  for  the  base  process  to  be  unstable  with  respect  to  all  but  a  very  small 
subset  of  the  initial  disturbances  under  consideration. 

Now,  to  show  that  u,  <  0  is  a  necessary  condition  for  the  base  state  to  be 
linezirly  imstable  with  respect  to  any  perturbation  in  the  full  class  introduced  in 
Section  3.2  consider  the  dependence  of  the  following  energy  S  on  time: 


-t-oo  0 


^(0  =Y  J  J  +  b]v,l{xi,X2,t))  dxydX'i 


+  00  -foo 


+Y  J  J  {v^{xuX2,t)  +  {al-vl)v,l{xi,X2,t)  +  blv,l{xi,X2^t))dxidx2 


—  oo  0 


V<€  [0,<.). 


(3.8.5) 


In  (3.8.5)  ^  is  a  positive  constant  which  ceirries  imits  of  length  and  the  function 
V  :  ]R^  X  M+  —*  Mis  defined  via 

w{xi,X2^t)  =  v{xi  —  vot,X2^t)  V(xi,i2,<)  E  X  jR+.  (3.8.6) 


By  (3. 5. 5)], 2  and  the  stipulated  square  integrability  of  and  it  is  evident 
that  S{0)  exists.  Assume,  then,  that  there  exists  a  positive  time  t,,  which  may 
be  very  small,  such  that  the  integrals  which  define  £  exist  on  the  time  interal 
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(0,<*).  Recalling  (3.1.5),  (3.3.8)  and  (3.3.10),  it  is  clear  that  S  is  non- negative  on 
its  domain  of  definition.  Now,  an  alternate  definition  of  linear  stability  in  terms 
of  the  energy  S  is  that  it  remain  bounded  on  JR-^.  A  series  of  calculations  which 
use  (3.5.2)  and  (3.5.4)-(3.5.6)  then  show  that  the  ■power  t  is  given  by 

-J-OO 

£(t)  =  —piv^  J  s^{x2^^)  dx2  V<  €  [0,  <*).  (3.8.7) 

—  OO 

Certainly,  if  u*  >  0  then  £{t)  <  0  for  all  t  in  [0,  t*);  furthermore,  under  these 
circumstances,  the  interval  over  which  S  is  defined  can  be  extended  in  increments 
to  leading  to  the  following  inequality 

£{t)  <0  V<  €  m+.  (3.8.8) 

The  condition  u*  >  0  is,  therefore,  sufficient  to  ensure  that  the  energy  S  re¬ 
mains  bounded  for  all  time  and,  hence — according  to  the  above  definition  of 
linear  stability — that  the  base  process  is  stable  with  respect  to  all  perturbations 
of  the  type  introduced  in  Section  3.2.  Now,  since  v,  ^  0  by  assumption,  it  is  clear 
that  u*  <  0  is  a  necessary  condition  for  the  beise  process  to  become  unstable  with 
respect  to  etny  perturbation  of  the  type  under  consideration. 

From  the  foregoing  discussion  it  is  cleeir  that  admissibility,  the  assumed 
smoothness  of  V  or  and  (3.1.9)  imply  that  when  uq  =  0>  the  base  driving 
traction  /o  =  0  and  V'(0)  >  0  or  ^'(0)  >  0  and,  hence,  a  mechanically  equili¬ 
brated  base  process  of  the  kind  introduced  in  Section  3.1  is — zis  in  the  inertia-free 
case  st  wxtli  ct  to  all  perturbations  considered  in  this  work. 

To  recapitulate,  the  calculations  performed  this  section  show  that  u*  <  0  is  a 
necessary  condition  for  the  base  process  to  be  unstable  with  respect  to  any  initial 
disturbance  of  the  type  under  consideration  and  a  sufficient  condition  for  the 
base  process  to  be  unstable  to  all  but  a  small  subset  of  these  initial  disturbances. 
Comparing  the  results  obtained  in  Section  3.6  with  those  obtained  here,  it  is 
apparent  that  inertizd  effects  do  not  significantly  effect  the  linear  stability  criteria. 
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3.9.  Discussion.  The  analysis  of  Sections  3.6  and  3.8  shows  that — 
regardless  of  whether  intertial  effects  are  included  or  not — u,  <  0  is  necessary 
for  the  base  process  to  be  linearly  unstable  with  respect  to  any  perturbation  of 
the  type  purveyed  in  Section  3.2  and  sufficient  for  the  base  process  to  be  vuistable 
with  respect  to  all  but  a  small  subset  of  these  perturbations.  In  Section  3.7  it  is 
shown,  in  the  absence  of  inertial  effects,  that  if  u,  <  0  then  the  linear  instabihty 
will  manifest  itself  in  a  manner  whereby  the  morphology  of  the  interface  evolves 
so  as  to  develop  plate-like  or  dendritic  structures.  Recall,  from  the  alternate  def¬ 
initions  (3.4.13)  and  (3.4.14)  of  u,,  that  v,  <  0  can  occur  only  if  the  relevant 
kinetic  response  function  V  or  9  is  locally  decreasing  at  /o  or  vq,  respectively. 
Is  it  physically  plausible  for  F  or  to  display  such  non-monotonicity?  Recall 
from  Section  2.3  that  admissibility — from  the  perspective  of  the  Clausius-Duhem 
inequality — does  not  restrict  the  monotonicity  of  the  kinetic  response  function. 
Owen,  Schoen  &  Srinivasan  [26]  suggest,  moreover,  that  unstable  kinetics  of 
the  sort  where  <p  has  a  single  maximum  as  a  function  of  Vi, — and  thus  must  be  a 
non-monotonic  function  of  Vn — may  be  responsible  for  the  rapid  growth  of  plate¬ 
like  structures  which  is  observed  experimentally.  Furthermore,  there  exist  other 
physical  contexts,  the  most  notable  of  which  include  imstable  crack  growth  and 
the  slip-stick  peeling  of  tape,®  where  the  analogues  of  such  non-monotonic  kinetic 
response  functions  are  considered  physically  acceptable. 

It  is  reasonable,  based  on  the  foregoing  discussion,  to  refer  to  the  type  of  linear 
instability  which  occurs  when  u*  <  0  as  kinetic  instablity.  This  investigation  has 
demonstrated  that,  under  the  current  kinematical  and  constitutive  restrictions,  in 
a  purely  mechanical  context,  independent  of  whether  inertial  effects  are  accounted 
for  or  not,  the  only  means  by  which  a  linear  instability  involving  the  emergence  of 
plate-like  or  dendritic  structures  from  a  planar  interface  can  occur  is  if  a  kinetic 
instability  is  present.  It  is  possible  that,  in  a  broader  context,  other  brands  of 
instability  may  be  present.  This  may,  in  particuleu:,  be  true  when  thermal  effects 

®  See  Aifantis  [7],  Aubrey,  Welding  &  Wong  [8]  and  Maugis  &  Barquins  [23]. 
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axe  taken  into  consideration.  An  investigation  which  takes  both  mechanical  and 
thermal  effects  into  consideration  is  performed  by  Fried  [15]. 
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APPENDIX 

In  this  appendix  (3.7.1)  and  (3.7.3)  are  established.  Consider  the  inertia-free 
initial  value  problem  comprised  of  (3.5.7)-(3.5.9),  (3. 5.2)2,  (3. 5.5)3  and  (3.5.6). 
Define  v  :  IR^  x  1FI+  IRhy 

w{xi,X2,t)  =  v{rji{xi,t),r}2(x2),t)  \/{xi,X2,t)  e  x  1R+,  (A.l) 


where  J7i(-,  t)  :  IR  JR  2md  7)2  IR  —*  IR  are  given  by 


rji{xi,t)  = 


{xi  —  vot)  if  xi  <  vot, 
(xi  —  Vot)  if  xi  >  Vot, 


(A.2) 


for  each  t  in  M+,  and 


Tj(x2)  =  X2  Vx2  €  ]R. 


(A.3) 


Then,  in  terms  of  v,  (3.5.7)-(3.5.9)  yield,  for  each  t  in  1R+, 

u,oa  =  0  on  ((— oo,0)  U  (0,  00))  X  iR, 
Or^rU,!  (0-|-,-,t)  —  a/6(U,i  (0— ,-,t)  =  0  on  ]R, 

v{0+,-,t)-v{0-,-,t)  =  {‘yi-yr)s{-,t)  on  M 


(A.4) 


OrbrV,!  (O-F,  •,  t)  +  a/fe/u,i  (0-,  •,  t)  ~ 


Now,  from  (A.4)2,4  it  is  clear  that 


2t;* 


(7/  -  7r) 


s{-,t)  on  IR. 


u,!  (0-f,  •,#)  -  y,i  (0-,-,0  = 


v^,{aibi  -  arbr) 
aibiarbriji  -  7r) 


i(-,f)  on  IR, 


(A.5) 


for  each  t  in  1R+  Consider,  next,  (A.4)i,3  and  (A.5).  Since  v(-,-,t)  is  harmonic 
on  ((-00, 0)  U  (0,oo))  X  JR  for  each  t  in  M+,  the  jumps  in  the  normal  derivative 
of  V  and  in  v  itself  across  the  line  I  —  {(»7i,»72)|  Vi  —  ^•>V2  ^  -K}  indicate  that  it 
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can  be  represented  by  the  sum  of  appropriate  single-  and  double-layer  potentials 
with  densities  proportional  to  i  and  s,  respectively.  That  is,  for  each  t  in  JR^, 


+  00 

v{rii,rj2,t)  =  J  iLi(T]i,T)2-C)HC,i)  +  L2(vi,‘n2-C)siC,'t)) 

— oo 


V(T/i,r?2)  e  ((-oo,0)  U  (0,oo))  X  iR,  (A. 6) 


where  La  :  \  {(0, 0)}  x  iR+  — »  ]R  are  given  by 

Liini.m)  =  ^  7r)  elR^\  {(0,0)}, 

L^(nur)l)  =  ^3^  €  B?  \  {(0,0)), 


(A.7) 


and  where  is  as  defined  in  (3.6.6).  With  the  aid  of  (A.1)-(A.3),  (A. 6)  gives 
the  representation  (3.7.1)  of  xa.  Next,  use  of  standard  identities  from  potential 
theory  gives  the  following  expression  for  the  limiting  values  of  u,i  (•,  •,<)  on  either 
side  of  I 


V,l  (0±,T72,f)  = 


-f-OO 

1  fs,2{C,i)dC  V^jaibl  -  Orbr)  .. 

2x7  J72  -  c  2aj6jar6r(7i  -  7r)^^ 

— oo 


^t]2  €  ]R,  (A. 8) 


for  each  t  in  1R+.  Substitution  of  (A.8)  into  either  (A. 4)2  or  (A.4)4  then  yields  the 
fimctioned  equation  (3.7.3)i.  The  initial  condition  (3. 7. 3)2  follows  directly  from 
(3.5.5.)3. 
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Figure  1:  Graph  of  the  shear  stress  response  function  r. 
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Figiire  2:  Graph  of  the  kinetic  response  function  V. 
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Figure  4:  Superimposed  graphs  of  the  h  and  «(•,<')  for  some  <'  ic- 
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